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We investigate cosmological particle production in spacetimes where Lorentz invariance emerges in the in- 
frared limit, but is explicitly broken in the ultraviolet regime. Thus these models are similar to many (but not all) 
models of quantum gravity, where a breakdown of Lorentz invariance is expected for ultraviolet physics around 
the Planck/ string scale. Our specific model focuses on the boost subgroup that supports CPT invariance and 
results in a momentum-dependent dispersion relation. Motivated by previous studies on spacetimes emerging 
from a microscopic substrate, we show how these modifications naturally lead to momentum-dependent rainbow 
metrics. 

Firstly, we investigate the possibility of reproducing cosmological particle production in spacetimes emerging 
from real Bose gases. Several papers have been written on the analogy between the kinematics of linearized 
perturbations in Bose-Einstein condensates and effective curved-spacetime quantum field theory. Recently, in 
[Phys. Rev. A 76 (2007) 043625] we have studied the influence of non-perturbative ultraviolet corrections in 
time-dependent analogue spacetimes, leading to momentum-dependent emergent rainbow spacetimes. We show 
that models involving a time-dependent microscopic interaction are suitable for mimicking quantum effects in 
FRW spacetimes. Within certain limits the analogy is sufficiently good to simulate relativistic quantum field 
theory in time-dependent classical backgrounds, and the quantum effects are approximately robust against the 
model-dependent modifications. 

Secondly, we analyze how significantly the particle production process deviates from the common picture. 
While very low-energy modes do not see the difference at all, some modes "re-enter the Hubble horizon" during 
the inflationary epoch, and extreme ultraviolet modes are completely insensitive to the expansion. The analysis 
outlined here, because it is nonperturbative in the rainbow metric, exhibits features that cannot be extracted 
simply from the standard perturbative modification of particle dispersion relations. However, we also show 
how the final result, after many e-foldings, will approach a time-independent exponentially decaying particle 
spectrum. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



and 



The use of Bose-Einstein condensates (BECs) as analogue 
models for quantum field theory in spatially flat k — Q 
Friedmann-Robertson- Walker (FRW) geometries 



ds^ = gab dx°- dx = -dr^ 



1=1 



(1) 



(in d spatial dimensions) has recently been extensively ex- 
plored in [IJ. The proposal GISIllEllSlTllSlis based on 
the analogy between the equation of motion for collective ex- 
citations around a macroscopically occupied ground state of 
an ultra-cold weakly interacting gas of Bosons {i.e., the novel 
state of matter referred to as the Bose-Einstein condensate), 
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(2) 



the covariant free-field equation (Klein-Gordon equation) for 
spin-0 massless particles. (The indices a, h run from to 
d whereas i, j run from 1 to d for the spatial coordinates 
only). Here is a symmetric covariant rank two tensor 
whose entries include purely collective (mean-field) variables 
c = c(t,x), the speed of sound, and v = v(t,x), the back- 
ground velocity. Using g"^^ g^c — <5"c we are able to define an 
effective line-element for the mean-field, 

ds^ = (-) ^ [-(c2 - dt" - dt^^ dx + dx^ dx] . 

(3) 

Here the quantity U arises from the microscopic description of 
a Bose-Einstein condensate, and represents the inter-atomic 
potential. 

It has previously been shown, see [TOl, that Eq. (|2]i is 
an equation for perturbations in the collective variables, the 
phase 6 and its conjugate momentum 11^ (see also 12] |3] 
HI El m El HI). This kind of behavior is not restricted to 
Bose-Einstein condensates. Indeed the first modern paper 
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on analogue models for gravity focussed on ordinary fluid 
mechanics; see 18j. Since then numerous media have been 
suggested as substrates to develop analogue models for grav- 
ity 1 11] [12] [m T?!. Nevertheless it is sometimes difficult to 
isolate the fundamental principles behind the analogy. 

Inspired by the Bose gas, where both the microscopic and 
macroscopic theory is understood to an adequate extent, we 
suggest that it might be useful to transfer the experience 
gained from the various known examples for emergent space- 
times into a more general framework: 

Emergent spacetimes involve (i) a microscopic 
system of fundamental objects {e.g., strings, 
atoms, or molecules); (ii) a dominant mean- 
field regime, where the microscopic degrees 
of freedom give way to collective variables; 
(iii) a "geometrical object" {e.g., a symmetric 
tensor) dominating the evolution for infrared 
classical and quantum excitations around the 
mean-field. 

Within certain limits we are free to choose our geometrical 
object to mimic Einstein's theory of gravity, {e.g., a symmet- 
ric rank two tensor gab that is conformal to an exact solution 
obtained by solving the Einstein equation for a physically rea- 
sonable stress-energy tensor). 

Analogue models for gravity are emergent space- 
times that are specifically adjusted to mimic 
as closely as possible Einstein's geometrical 
theory of gravity. 

The fundamental difference between emergent spacetimes and 
general relativity becomes obvious if we consider the follow- 
ing statement: Emergent spacetimes, as they appear in the 
analogue model programme, are "a short and simple way" to 
summarize the kinematics of linearized perturbations in a geo- 
metrical sense, without requiring the notion of a stress-energy 
tensor or the Einstein equations. 

Thus, the analogy is (currently) restricted to the kinematic 
behavior of the system, and any back-reaction between the 
excitations and the mean-field might not be in analogy with 
Einstein's theory of gravity. Therefore it is advisable to re- 
strict the use of analogue models for gravity as toy models for 
semi-classical quantum gravity at the level of curved space- 
time quantum field theory, where the gravitational field is a 
purely classical field. 

Indeed, analogue models for gravity are very idealized 
constructions, and require careful treatment with respect to 
the model-dependence to obtain the desired curved-spacetime 
quantum-field-theory effects. This has first been pointed 
out in ifTSl [T6l . where the authors analyzed the robustness 
of Hawking radiation against modifications in the excitation 
spectrum (most relevant to fluid dynamics: super- and sub- 
luminal dispersion). 

In a related previous work fT] we made use of classical 
phase space methods to numerically study "cosmological par- 
ticle production" in a realistic BEC. In principle, there are two 



different ways to mimic an expanding universe in a Bose gas. 
In il] B l the authors studied a freely expanding condensate 
cloud, such that the density is a function of time. This idea has 
been further investigated in |7|. In both — our previous pa- 
per HI and in this work we pursue an alternative (more recent) 
idea introduced in ||2][3l, where the authors achieved FRW ge- 
ometry through a time-dependence in the effective two-body 
inter-atomic potential. (A change in the inter-atomic poten- 
tial is related to a change in the speed of sound, and con- 
sequently corresponds to a time-dependent acoustic metric.) 
As expected our simulations verified the analogy (in the hy- 
drodynamic limit) to a good extent. However, it is crucial 
to include quantum pressure effects to understand the ultravi- 
olet /trans-phononic part of the quasi-particle spectrum|65|. 
Thus we conclude that our specific model — an ultra-cold gas 
of Bosons with time-varying atomic interactions — shows a 
specific model-dependence that affects the high-momentum 
particle production process. One lesson that can be drawn 
from this new insight is that a careful choice of condensate 
parameters is necessary in order to mimic the desired curved- 
spacetime quantum field theory effect. 

The main focus in this article is to understand the differ- 
ences between ordinary quantum field theory in a "real" clas- 
sical FRW-background and quantum field theory in an emer- 
gent FRW-background in a Bose gas. This question has been 
motivated by many effective field theories, where Lorentz in- 
variance is broken at ultraviolet energies in a similar manner 
as in the emergent spacetimes we are investigating. 

Emergent spacetimes exhibit an emer- 
gent / ejfective Lorentz symmetry for 
low-energy / infrared excitations around 
the macroscopic field. This symmetry will be 
broken in the high-energy / ultraviolet regime, 
that is at scales dominated by the underlying 
microscopic theory. These corrections are of 
a non-perturbative nature. 

For now, we leave general questions aside, and focus on 
our specific emergent spacetime model, investigating the in- 
fluence on cosmological particle production due to ultraviolet 
corrections. We include quantum pressure effects — which 
lead to a nonlinear excitation spectrum — as nonperturba- 
tive ultraviolet corrections to the emergent geometry. We will 
show that this leads to momentum-dependent "rainbow space- 
times". Partly motivated by extant literature |17, 18 1 and our 
numerical results from the simulations of quantum effects in 
realistic Bose-Einstein condensates, we present a coherent in- 
terpretation for the quasi-particle spectrum produced in FRW 
rainbow spacetimes llT8l[T9ll20l 12111221123 1. 

The "FRW rainbow metrics" recover the "standard" FRW 
type geometries in the infrared (phononic) regime, but ex- 
hibit momentum-dependent modifications in the ultravio- 
let (trans-phononic) regime. We show that in such a ge- 
ometry the scale factor for the universe a{t) is effectively 
momentum-dependent a{t) ak{t), consequently leading to 
a momentum-dependent Hubble parameter: H{t) Hk{t). 
In addition, we are dealing with emergent spacetime that ex- 
hibits a time-dependent effective "Planck-length". Ultimately, 
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the particle production process in our specific analogue model 
will show deviations from the one expected in "standard" 
curved-spacetime quantum field theory. 

For the "analogue gravity programme" it is important to 
check the robustness of cosmological particle production 
against those modifications. Is there a window where the 
analogy is good enough, for example, to use the Bose gas as 
a toy model for inflation? 



quires a more accurate analysis at the level of the emergent 
geometry, leading to the concept of a rainbow metric. In other 
words we extend the treatment for ultraviolet modes from flat 
to curved spacetimes. 



Ultraviolet non-perturbative corrections: 
Quantum pressure effects 



In this paper we will show that in the present model the 
particle production process is in general not robust against 
the model-specific modifications. However, for short time 
duration expansion scenarios, inflation can be simulated in a 
realistic Bose-Einstein condensate. 

We perform a quantitative analysis to describe short time 
as well as long-lasting expansions. We use the ratio between 
the mode frequency and the Hubble frequency (the inverse 
of the rate of change in the size of the emergent universe), 
to show that our specific model indeed shows significant 
deviations in the ultraviolet regime, e.g., both crossing and 
re-entering of the "Hubble horizon" during the inflationary 
epoch; a{t) ~ exp(i). 

By means of our numerical analysis we are able to "read- 
off" the final spectrum after a sufficiently long-lasting expan- 
sion. We explain why the final particle spectrum — after 
an infinitely long-lasting inflationary epoch — will be time- 
independent and finite. 



II. EMERGENT SPACETIMES: 
EXCITATIONS IN BOSE-EINSTEIN CONDENSATES 

In the following we extend the calculations presented in our 
previous work 1 1 1. The intent is to derive an emergent geome- 
try that is able to simultaneously represent both the phononic 
and trans-phononic excitations in a Bose-Einstein condensate. 
It is a well established result — both in theory and experi- 
ment — that the energy-momentum relation for excitations in 
Bose-Einstein condensates is given by the nonlinear Bogoli- 
ubov dispersion relation. In|18 19, 20, 21 , 22, 23 1 it has been 
shown, (within the analogue model for gravity point of view), 
that phononic modes are relativistic modes, since they exhibit 
Lorentz symmetry. The corresponding emergent geometry is a 
flat spacetime with Lorentzian signature (Minkowski metric). 
Higher energy (trans-phononic) modes break Lorentz symme- 
try, the excitations being "supersonic". The full dispersion re- 
lation for excitations in a realistic Bose-Einstein condensate 
is similar to that expected to occur in certain effective field 
theories. It is also known that this kind of modification is non- 
perturbative, since it originates in density fluctuations visible 
only at small scales (in the order of the healing /coherence 
length). These modes are trans-phononic modes, where quan- 
tum pressure effects are no longer negligible. We show that 
in a Bose-Einstein condensate with time-dependent conden- 
sate parameters, the crossover between phononic and trans- 
phononic modes is also time-dependent, and therefore re- 



in IT] we derived a coupled pair of equations. 



dth + V ■ 



\ m 



V0 +(nov) 



0, 



U 



dte + v -^9+ =r h = 0: 
n 



(4) 
(5) 



for quantum fluctuations in the condensate density fi and 
phase 6. The background velocity v is given by 



V = - V^o , 
m 



(6) 



as the gradient of the condensate phase Oq. Here m is the mass 
of the fundamental Bosons, uq the number density, and h the 
Planck constant. The differential operator U has been defined 
asdllll 



where the differential operator 



D2 - 



2 \ nl nl 



no 



(7) 



(8) 



accounts for the first-order correction obtained from lineariz- 
ing the quantum potential term whereas the inter-atomic po- 
tential is explicitly given by 



U 



(9) 



The scattering length ascatt(i) represents the s-wave scat- 
tering term. (See, for instance \XL ,24], and the discussion 
above.) Note the tilde notation is used to emphasize the fact 
that we are dealing with a differential operator 

To develop the analogy between condensed matter physics 
and curved-spacetime quantum field theory, see Eq. (j2]i, it is 
necessary to combine Eq. (|4]) and Eq. (|5]l into a single equation 
for the phase fluctuations. For that it is important to note that 
we can rearrange (|5]l to make h the subject. That is 



= hu- 



■ V6I 



= hu 



_i D9 



(10) 



where Dd/Dt is the rate of change of 6 following a small 
volume of the fluid, the material derivative of 0. If the fluid 
is at rest, the material derivative reduces to D9/Dt — ^ dtO, 
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but for a moving fluid it represents a "fluid-following" deriva- 
tive. The integral differential operator U^^ cw formally be 
expanded as 



2m 



(11) 



2m 



n 3 



2m 



U-^D2 U-'^D2 U-^D2 U-^ 



where the formal series converges only on the subspace of 
functions spanned by the eigenfunctions whose eigenvalues 
satisfy 



A' 



2to 



U-'D2 < 1 



(12) 



Since D2 and U are second-order linear differential operators, 
the inverse always exists as an integral operator (that is, 
in the sense of being a Green function). Expanding this Green 
function as in Eq. 
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above is a convenience that allows us 
to interpret U^^ as a sum of differential operators, but this is 
not a fundamental limitation on the formalism. 



1. Acoustic and rainbow metrics 



We use equation (10 1 to eliminate h completely in our equa- 
tions of motions, by substituting it into equation Q. That 
yields a single equation for the perturbed phase 



(13) 



where we have introduced the matrix /"'' with inverse- 
differential-operator- valued entries: 



rab 



(14) 



Note that in general U^^ is an integral operator so Eq. ( 13 1 
is an integro-differential equation. If we additionally require 
that there exists an (inverse) metric tensor gab such that 



(15) 



where g is the determinant of gab, then the connection is 
formally made to the field equation for a minimally coupled 
massless scalar field in a curved spacetime; see Eq. 
Now in many situations of physical interest, the differential 
operator U can be usefully approximated by a function — for 
instance the hydrodynamic and eikonal limits. 



a. Hydrodynamic approximation: In the hydrodynamic 
limit where the quantum pressure is neglected one has 



so that 



\Uno\ > \{hy2m) D2 no\ 



U ^U. 



(16) 



(17) 



Within this approximation we obtain 



rab ^ ^ 

^ u 



- 1 













(18) 



To mimic curved spacetime classical and quantum field the- 
ory effects the acoustic metric has to be formally in agreement 
with Einstein's theory of gravity. For this low-momentum ap- 
proximation we are able to define a common (for all wave- 
lengths k) speed of sound in the condensate; 



(19) 



Applying Eq. ( 15 1 we obtain. 



9ab 



cm 



_ noU{t) 




m 













(20) 



as the acoustic metric or analogue metric. In this case all 
collective excitations behave as sound waves with the usual 
linear dispersion form 



ujq ~ ck 

and the quanta of excitations are thus phonons. 



(21) 



b. Eikonal approximation: Let us now invoke a differ- 
ent approximation that holds for trans-phononic modes, so 
that the elements of the matrix f°-^ can be treated as (possi- 
bly momentum dependent) functions, rather than differential 
operators: Consider the eikonal limit where U can usefully be 
approximated by a function 



C/^[/fc(i,x)-[/(t,x) + 



Amnn 



(22) 



which we shall conveniently abbreviate by writing Uk- Be- 
yond the hydrodynamic limit we obtain 



h 


- 1 






- v' 


noUk Xii yiyi 




m 



(23) 



Note that in the eikonal approximation the k dependence hid- 
ing in Uk will make this a momentum-dependent metric, a so- 
called rainbow metric. The metric tensor is explicitly given 
by 



gab 



riQ h 
Ck m 



(24) 



5 



where we have introduced the quantity 



(25) 



for a uniform condensate density no(t, x) — tzq. It is conve- 
nient to define 



7qp = 



2m 



(26) 



where 7qp <C |c(i)/fc| is a useful indication for the hydrody- 
namic limit; compare with Eq. ([T6]l. The dispersion relation 
in the eikonal limit is 



(27) 



and hence violates "acoustic Lorentz invariance" . This is not 
surprising at all, since we know the quasi-particles become 
"atom-hke", and so non-relativistic, at high momentum ||23]| . 

c. Acoustic or rainbow spacetimes? It is only in the 
acoustic/hydrodynamic limit that the k dependence of Uk and 
Cfc vanish, so that the rainbow metric is reduced to an ordinary 
Lorentzian metric. 

An interesting consequence of the Bogoliubov theory in 
Bose condensates is that in general the excitation spectrum 



displays nonlinear dispersion (see Eq. (27 1), being linear {i.e., 
phononic) for low |k| and becoming quadratic {i.e., free- 
particle like) at large |k|. When the nonlinear dispersion (27 1 



is incorporated into analogue models of gravity it is equivalent 
to breaking Lorentz invariance fTTI [191 ifol 1211 122ll24ll . 

The hydrodynamic approximation is a statement about the 
smallness of the quartic term in the dispersion relation for 



Lof. (|27|), with respect to the quadratic term. The usual line of 



argument is that the smallness of 7qp makes it possible to ne- 
glect the second order in Eq. ( 27 1 for low -energy excitations, 
where 7qp|fc| <C c{t). Here Lorentz invariance is an emergent 
symmetry. However, it is important to realise that the propa- 
gation speed can be a function of time, c — c{t), and hence 
if c(<) then one is dealing with a system that eventually 
violates Lorentz invariance at all energy scales. In principle, 
there are no theoretical or experimental restrictions to prevent 
U{t) (X c{t) becoming arbitrarily small. 

In the specific cases we are interested in, we are confronted 
with exactly this situation, and therefore a more subtle 



analysis is required as to whether the acoustic metric ( 20 1 is a 
sufficient approximation, or whether we have to use the more 
sophisticated concept of a rainbow spacetime (24i. We will 



elaborate on the point in Sec. IV 



d. Limitations of the rainbow analogy? It is worth re- 
iterating the assumptions used in making the above anal- 
ogy. The theory resulting from linearized quantum fluctua- 
tions leads to a free field theory — that is, the modes of the 
field 9 are non-interacting. Since we are using a linearized 
theory, the interactions between quantum fluctuations them- 
selves, and between quantum fluctuations and the condensate 
mode are neglected. This is equivalent to assuming the metric 
tensor is an externally-specified classical quantity. 



Our numerical simulations carried out in [Tl do not require 
these assumptions; there all modes of the system are included 
and these modes are able to interact — albeit weakly — via 
the nonlinear interaction term. Thus we were able to explore 
the validity of the assumptions of the free-field theory Iii25ll . 

We further note that while the present form of the analogy 
only holds for massless scalar (spin zero) particles, in gen- 
eral it is possible to modify the formalism to include massive 
minimally coupled scalar fields at the expense of dealing with 
more complex EEC configurations, e.g., a two-component 
EEC ||T8l|T8l[l3|20l|2l]|^|26l|27|. In EEC language one 
would explain this situation in terms of a dispersion relation 
with a gap. 



2. Commutation relations 

To derive the analogy presented above, we approximated 
and transformed our field operators several times, see ||9| 
and ID . These are canonical transformations preserving com- 
mutation relations: 

1. ■0(t,x) and'(/>t(t,x): The single Eoson annihilation and 
creation operators; where 

[7/;(t,x),7/;(i,x')] = 0, (28) 
[7/;t(t,x),V^t(t,x')] = 0, (29) 
[V;t(t,x),vJt(t,x')] = <5(x-x'); (30) 

2. 5il}{t,yL) and 5ijj^{t,yL): Decomposition into a single 
coherent mode ip{t,x.) — (?/'(<, x)), and the quantum 
excitations Sil;{t, x) around it. Altogether, ^{t, x) = 
?/;(i,x) + (5V;(t,x) andV;t(t,x) = ij*{t,x)+Sjp'<{t,x), 
where 

[Sij{t,x),6i;it,x')]^0, (31) 
[57/;t(t,x),5V^t(t,x')] = (32) 
[,57/;t(t,x),5V^t(t,x')]-<5(x-x'); (33) 

3. n and 6: Mapping onto Hermitian phase and density 
fluctuation operators, as studied (for example) in |2 3|. 
Here we made use of the fact that the macroscopic field 
V^(x) is complex and so for topologically trivial regions 
— without zeros or singularities — one can always ex- 
press it as ijj{t, x) = y/n{t, x) exp{i9{t, x)). Lineariz- 
ing around the two parameters of the complex-valued 
field, 9 ^ 9o + 9 and n ^ uq + n, we can write 
rp ip + y/riQ ^ 2^ + 'i9j, and its Hermitian conju- 
gate, such that 



[n(i,x),n(t,x')] -0, (34) 
^(t,x),^(t,x')j =0, (35) 

h{t, x), ^(i, x')j =iS(x- x') . (36) 
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4. Finally, we are able to use the equation of motion ( 10 1 
to formally express h in terms of 9: 



1 De{t,yi) 1 De{t,y.') 
U Dt 'u Dt 



1 De{t,y.) 



U Dt 



= 0, 
,0(t,x') 



= 0, (37) 
(38) 

I 5(x - x') ; (39) 



There are two different ways to view the resulting commutator 
relations in the hydrodynamic and eikonal limit: A condensed 
matter point of view in terms of time-dependent commuta- 
tors, or a commutator relationship for the phase perturbation 
6 and its conjugate momentum 11^ on the emergent spacetime. 

a. Condensed matter point of view: First we face the 
problem of how to deal with the differential operator D2, 
which involves studying two interesting limits where the 
commutation relation takes a simpler form. 



In the hydrodynamic approximation we get 



D6'(i,x) 

Dt ' 



6{t,^') 



m 

ih 



5(x-x'), 



(40) 



which is now a time-dependent commutation relation. 



For U{t) the hydrodynamic commutator vanishes 
completely and we are left with purely classical statements 
for 9. This situation changes significantly if one instead 
considers the eikonal approximation. 

In the eikonal approximation we get (in momentum space) 



D9{t,k) 
Dt ' 



(i,k') 



U{t)- 



: Ofefe' , (4i) 



and the commutator does not vanish for U (<) 0, though it 
can vanish if U {t) becomes negative. 

This suggests that the presence of D2 cannot in general be 
neglected for a time-dependent atomic interaction U{t). 

b. Emergent spacetime point of view: An alternative in- 
sight can be gained if we define an emergent Lagrange density. 



1 



f'''da9db9, 



(42) 



in correspondence with Eq. The momentum conjugate to 
9 is given by 



n, := 



dC 

d{dj) 



= -f'db9, 



(43) 



and hence we evaluate the conjugate momentum to 9 as, 

h D9 

- ^ — . (44) 

U Dt 



With this new insight we are able to add another set of com- 
mutation relations, one that makes only sense after having in- 
troduced the emergent spacetime: 

5. The phase and density operators are a canonical set of 
quantum field operators and conjugate field operators, 



0(i,x),^(i,x') 



0, 

= 0, 



(45) 
(46) 



n,^(i,x),n^(i,x') 

0(t,x),n,^(i,x')] =^5(x-x'); (47) 
in an effective curved spacetime represented by 



Eq. (20 1, for a massless spin-zero scalar field. 



In the hydrodynamic limit U ^ U we recover 

h D9 



the standard result for the conjugate momentum in curved 
spacetime; for more details see Il28ll29]| . 

The conjugate momentum /or ultraviolet modes U ^ Uk 
is given by 



h D9 
lADt 



(49) 



as the conjugate momentum in our rainbow geometry. 



This is a significant result, since it shows that knowledge of 
the emergent spacetime picture provides a deeper insight into 
the full dynamics for the density and phase perturbations, and 
explains the explicit time-dependence in their commutation 
relations. 



B. FRW- rainbow geometries: 
Specific time-dependence for atom-atom scattering 

Clearly time dependence can enter in any of the parameters 
no, c, and v. We now focus on the case where v = 0, and no 
is constant throughout space and time, so that the system is 
homogeneous. This choice of parameters leads to the specific 
class of A: = spatially flat FRW spacetimes ll3l l30ll . Such 
geometries are always conformally flat and at any particular 
time the spatial geometry is simply that of flat Euclidean 
space. All the time dependence is contained entirely in 
the speed of sound given by Eq. 



19 1. In a homogeneous 



condensate the speed of sound and the scale factor are 
position independent, and within the acoustic / eikonal limit 
a separation of the field operators into time and position 
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dependent parts is possible. 

We introduce the dimensionless scale function b{t) so that 
the interaction strength (or equivalently the scattering length) 
becomes time-dependent. 



U{t)^Uob{t), 
and under consideration of Eq. (|7| we get 

u = uom, 

and a time-dependent differential scaling operator 6(t); 



m = b{t) 



2m U, 



Do 



(50) 



(51) 



(52) 







It is important to realize that a change in the interaction 
strength, U U{t), inevitably involves a shift in the 
crossover between the phononic and trans-phononic regime, 
and hence the nature of collective excitations in the conden- 
sate. To see this relationship more clearly we write down the 
eikonal approximation of Eq. ([52|, 



Uk{t) = Uo bk{t) = Uo {b{t) + {k/Kf) , (53) 
where for further convenience it is useful to introduce K, 



K 



7qp 



(54) 



which represents the crossover between the phononic and 
trans-phononic regimes. Here we defined cq = c(to), as 
the initial speed of sound. We would like to stress the im- 
portance of the last two equations for the understanding of 
everything that follows below. Motivated by the ultraviolet 
deviations in our results obtained by the numerical simula- 
tion of time-dependent spacetimes in a realistic Bose-Einstein 
condensate [l], we set out to find a description that includes 
trans-phononic modes into the spacetime picture. This forced 
us to generalize acoustic metrics to momentum-dependent 
(rainbow) metrics, that cover a larger fc-range. These rain- 



bow metrics, see Eq. ( 23 1, show a time- and momentum- 
dependent scaling in the effective scale function for the in- 
teraction strength. To show this we take &(io) = 1 at some 



arbitrary initial time <o- From the convergence constraint (12 1 



for a uniform condensate at rest at t = 0, and working in the 
eikonal approximation. 



Ifcl < \K\ 



(55) 



This implies that for modes k > K the spacetime picture be- 
gins to break down. This strongly suggests that one should 
consider 1/K as the analogue Planck length, ^pianck ^ l/K. 
For modes with wavelength k <^ K we recover the "stan- 
dard" geometry (i.e., momentum-independent spacetimes). 
These modes are phononic modes. Making use of our rain- 
bow geometries we extended the emergent spacetime pic- 
ture for higher energetic (ultraviolet) modes with wavelengths 
< k < K. For rainbow metrics the equation of motion 
for modes with wavelengths k > K sue most easily derived 



in momentum space. This can immediately be seen when 
one Fourier transforms the initial quantum fluctuations Stp and 
Si/j^, since then it is not necessary to introduce gab at all. The 
situation is more complicated if we involve time-dependent 
atomic interactions. At any later time t > Owe get a different 
convergence constraint. 



\k\ < 



(56) 



and with it a time-dependent limit on the breakdown of the 
validity of the emergent geometry. 



Planck 



(57) 



In other words, the analogue Planck volume ^pianck(i)'^ is 
— depending on the form of b{t) — shrinking or expanding. 
We will return to this point in Sees. IV C 2 and[V] when we 
explain particle production in a de Sitter-like universe, with 
a growing Planck volume. Note that for a time dependent 
condensate density the problem is different, and has been 
studied in El m 

It should be noted that in the acoustic approximation the 
scale factor. 



Uk{t) -> Uo b{t) 



(58) 



is indeed momentum-independent as expected. In practice 
a variation in the interaction strength is possible by using a 
Feshbach resonance Il2ll3l [3n[32]| . 



From Eq. ( 25 1 we obtain the time-dependent speed of sound 
in the eikonal approximation. 



Ck{i) = Co y/bit) + kyK^ , (59) 



and thus 



where 



LJkit) = Luo Vb{t) + kyx^ = LOO Vh{i), (60) 



bkit)^b{t) + k^/K\ (61) 



Their acoustic counterparts are obtained in the limit 
bk{t) b{t). 



The equation of motion ( 13 i, which depends on /" , there- 



fore shows no explicit dependence on the spatial dimensions 
d of the condensate. In contrast, the emergent spacetime 



-clbkitrdt^ + bkity 



(62) 



is explicitly dependent on d through gat, see Eqs. (20 1 



given as 



and p4| . Here the exponent a is dimension-dependent and 

d-2 

(63) 



d-1 
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In order to use a Bose-Einstein condensate (in the infrared 
limit) as an analogue model for Einstein's theory of gravity, 
we have to equate the Friedmann-Roberton-Walker universe 
line element given in Eq. ([T]i, with the line element ( 62 1. This 
explains why the choice of b{t) — to mimic a specific FRW- 
type universe scale factor a(r) — depends on the spatial di- 
mensions d. 



III. QUANTUM FIELD THEORY IN RAINBOW 
GEOMETRIES: 
PARAMETRIC EXCITATIONS 



Before we address problems of interest to the programme 
of cosmological particle production in a specific FRW-type 
universe, we prepare the necessary mathematical and physical 
formalism. Quantum field theory in rainbow geometries 
requires a careful treatment of the equation of motion and the 
commutation relations for the quantum field operators. There 
are three major differences between ordinary massless spin-0 
particles in real curved spacetimes and BEC quasi-particles in 
the particle production process. These are: (1) The numerical 
finiteness of the particle production in all circumstances 
due to the non-perturbative ultraviolet corrections, as we 
will show below. (2) We are taking the point of view that 
any Bose-Einstein condensate experiment will have a finite 
time duration, and therefore well-defined initial and final 
vacuum states. (3) Finally, there exists a preferred frame, 
the laboratory frame in which the experiment (or at this 
stage gedanken-experiment) is implemented. The laboratory 
time therefore can be viewed as the most relevant coordinate 
choice. 



Conditions (1) and (2) allow us to employ the method of 
instantaneous Hamiltonian diagonalization to adequately cal- 
culate the particle production. The condition (3), that our par- 
ticle detector is bound to the laboratory-frame simplifies the 
decision as to which time coordinate to choose. However, this 
limitation is also very disappointing. A major awareness re- 
sulting from "conventional" curved-spacetime quantum field 
theory has been the observer dependence of the particle spec- 
trum |33|: Different coordinatizations motivate different vac- 
uum choices. One might stretch the good nature of any con- 
densed matter experimentalist, asking for a co-moving par- 
ticle detector — adjusted to a suitable vacuum in the cho- 
sen coordinate system — in an infinitely long-lasting expand- 
ing emergent spacetime. Any unambiguous measurement re- 
quires one to stop the expansion, and to project onto a positive 
and negative plane wave basis for the final value of Uk{t). 



Linearized quantum excitations in condensate 

The equation of motion for quantum fluctuations in such a 
d-dimensional condensate are given by 

d, { J-dAt,^)] ^V20>,x) = 0; (64) 
\U{t) / m 



compare with Eq. ( |23| ). In the following we are going to 
rewrite the equation of motion in a more suitable form, ap- 
plying an auxiliary field, and simultaneously transform to mo- 
mentum space. 



1. Auxiliary field operators in Fourier space 

We use the differential operator R — U{t)/h to write the 
field operators in terms of auxiliary field operators x, where 
e = (^)i/2 ^. As long as [/ = U{t), and with it ^ = e{t) and 
X = is position independent, we can always write 




(65) 



The last statement is naively based on a Taylor series which 
only converges for modes with wavenumbers \k\ < \K\. 
Within this radius of convergence the transformation is exact, 
but can be extended to arbitrary /c-values in the eikonal ap- 
proximation, where U{t) — > Uk{t). The equation of motion 
for the mode operators Xk, 



V Uk{t) 

and the equal time commutation relations, 

[Xfc(t),Xfc'(t)] = 0, 
[dtXk{t).dtXk'{t)]^Q, 
[Xk{t),dtXk'{t)] = i5k,k' ; 



(66) 



(67) 
(68) 
(69) 



are now slightly more convenient than (45 1-(47 1. The function 
^k{t) is defined as 



nk{tf = clk' bk{t) 




l hk{t) 
2bk{t) ' 



(70) 



A. Particle production 



and therefore we have the connection between quantum field 
theory in FRW-type spacetimes and a parametrically excited 
harmonic oscillator. The notation ilfc(t) will become quite 
obvious for time-independent cases, where it reduces to the 
usual dispersion relation Vlk ~^ ujk', see Eq. (60 1. Notice the 
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overall factor, yJh/Uk{t), in Eq. (66 1. This is of importance 



for the particle production process in cases of discontinuous, 
and continuous (but not differentiable) changes in U{t). We 



will revisit this issue shortly in Sec. Ill A 5 a 



2. Mode expansion 



The equation of motion ( [66] l is a homogeneous differential 
equation, which can be written as 



where 



(71) 



(72) 



is a linear second order differential operator with a 2- 
dimensional solution space. 

A common tool in quantum mechanics is to describe quan- 
tum states — with unknown or variable number of particles 
— with respect to an orthonormal occupancy number basis 
{i.e., the Fock space basis in the infinite-dimensional function 
space, the Hilbert space of state). A Fock state is a quantum 
state consisting of an ensemble of excited non-interacting par- 
ticles. 



t \nk2 



I»kl,rik2,-) = , , , = («lcl)""(42) 

Vnki'"k2'--- L 

(73) 

with definite occupation numbers (riki, nk2, ■•■) in the modes 
(Xki,Xk2, •■•)■ The creation, d\, and annihilation, opera 



tors 



create or destroy a single-particle in the mode Xk; 



alii..., Tiki - 1, . 
Okil---, n^i, . 



= \/ri\d.\---,n\^i, ■■■) , 
= \/nki |---,n-kj - 1, •■•) 



(74) 
(75) 



The state |0) — short for |Oki , Ok2, •■•) — is a special state, the 
vacuum state. It is defined as the eigenstate of all annihilation 
operators with eigenvalue 0, such that ak\0) ~ 0. Thus 
any arbitrary quantum state \(fi) in the Fock spaces is a linear 
combination of all excited states. 



nk2, -- Fkl, «k2, •■ 



(76) 



that can be created out of the vacuum. Here Priki.nk?, ■■ 
the probability to measure the single Fock state with the mode 
occupation (nki,nk2, ••■)■ 

Within this framework we expand the mode operators in 
terms of destruction and creation operators. 



Xk{t) 



1 

71 



vl{t) dk + Vk{t) al^. 



(77) 



The coefficients Vk and are a set of linearly independent 
mode functions and any linear combination is a complete so- 
lution of 



Lvkit) = 0. 



(78) 



(Due to the isotropy of the s-wave scattering amplitude of 
the atomic interactions, the mode functions are also isotropic, 
i.e., Uk — Vk and uj^ — vl. It therefore seems unlikely 
that we would need to extend the specific analogue set-up to 
mimic anisotropic expansion scenarios. The case is different 
for analogue models involving a changing condensate density, 
no = nQ{t); where a non-uniform density expansion can eas- 
ily be achieved, for example see |4, 5 6 7|.) 

In order to obtain a "nice" canonical set of operators, 
which obey the common equal time commutation relations 



(seeEq. (45 i-(46 1), 



4 '4' 

afc,4' 



= 0, 
= 0, 



(79) 
(80) 
(81) 



the Wronskian W of the mode functions has to be normalized 
as follows. 



W[vk,vl] =vkvl- Vk v*k = 2i. 



(82) 



It is easy to see that — under the application of the equation 
of motion for the mode functions — the Wronskian is always 
time-independent. Note, a spin-zero scalar particle is its own 



We would also like to draw 



anti -particle, where {dk) 
attention to the fact that with a particular choice of Fock space, 
we select a particular subspace of the Hilbert space of states 
{i.e., the subspace that can be created out of the vacuum state, 
determined by our choice of dk). This will be of further in- 
terest when we discuss the validity and limitations of the Bo- 
goliubov transformation. 



3. Bogoliubov transformation 

In Eq. ( [77] i we can in principle pick some particular mode 
expansion, {vkY ■ Ak, with a specific set of mode func- 
tions vj^ = {vk{t),vl{t)), and set of mode operators A]^ — 
(flfc, a^j,). This choice is not unique and any other mode ex- 
pansion (li/i;)^^ • with mode functions = {uk{t),u*f^{t)), 
and mode operators = (bk, 6^^,) would have been possi- 
ble. The relation between these two different representations 
is referred to as Bogoliubov transformation. 

As mentioned above, the mode expansion {vkY ■ A cor- 
responds to an orthonormal basis in the infinite-dimensional 
function space, provided the Vk are normalized; that is 
W\vk{t),v1.{t)] = 2i for all times t. If we restrict our- 
selves to exclusively mapping between orthonormal frames 
— then the mode functions fulfill the normalization constraint 
W[uk{t), u*f,{t)] = 2i, and the mode operators obey the com- 
mutation relations given in Eqs. ( 79pT 



transformation is given hy Vk — M 
2x2 matrix, 

Pk ak 



Uk, 



— the Bogoliubov 
where M must be a 



(83) 
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The complex-valued coefficients Uk and are called Bogoli- 
ubov coefficients, and the transformation preserves the nor- 
malization condition if det(M) = 1, that is for 



\ak\ 



= 1. 



(84) 



At this stage we notice a minor technical subtlety when 
we look at summing over discreet k modes versus in- 
tegrating over continuous k modes. Equation (84i is 



appropriate for summing over discreet modes whereas 
for continuous modes we should strictly speaking use 
|afe|2 - |/3fe|2 = ^W(O) = V/{2ttY. Here ,5('^)(0) is the 
d-dimensional momentum space Dirac function. We will 
work with continuous integrals in momentum space but will 
suppress unnecessary occurrence of the volume V, where 
including it would lead to unnecessary clutter that might not 
aid understanding. 

The relationship between the "old" and "new" mode oper- 
ators, where (ufe)^ • A — {ukY ■ B and (ufe)^ = (wfc)^ ■ M\ is 
given as Ak = (A/t)-! • and Bu = (M^) • Ak- Here 

Mt = f "'^ and (Mt)-i ^( < \ , (85) 

\Pk al J \ -(3k ak J ' 

and therefore 

ak = albk - P*kb^_k, and a}_f^ ^ akb^_k - I3kh\; (86) 
bk = akcik + Pla}_f., and Slj. = a^al^, - /5feaj^ . (87) 

Obviously, the two vacuum states, ak |(a)0) — and 
ik l(b)O) = 0, are different, as bk |(a)0) for 13k ^ 0. 



4. Particles and lowest energy eigenstate 

Quasi-particles are excited states, and therefore they de- 
pend on the choice of (mode functions and the associated) vac- 



uum state. Under the application of Eq. (87 1 we can formally 
calculate the mean density of fe-particles, nT*), the occupation 
number for the mode Xfc^ with respect to the a- vacuum |(a)0) 
as follows: 



{ia)Q\blbk\ia)Q) 

Volume 



(88) 



The challenge is to find particular specified mode functions 
that represent physically meaningful vacuum states, that min- 
imize the expectation value of the Hamiltonian. Only then 
does the chosen vacuum correspond to the "actual" physi- 
cal vacuum, and its excitations describe "real" quasi-particles. 
(That is, we are looking for an unambiguous measure of par- 
ticles relative to some vacuum state.) The problem of find- 
ing the "right" vacuum state ("best" orthonormal frame) has 
been controversial from the outset. However, as we shall 
see, it is straightforward to find the zero-particle vacuum 
state for Minkowski spacetimes, and the whole controversy 



arises only for time-dependent cases, where the Hamiltonian 
H{t) — J-^^Ti-kit) is explicitly time-dependent; 



Hkit) = 



{[9t(^/C4(^)xfc(^))][a*(^/C4(^)x-fe(^))]} 



2Uk{t) 

■,^k{tfxk{t) X-k{t) ■ 



(89) 



In combination with the mode expansion ( 77 1 the Hamiltonian 
in momentum space is given by 

H{t) = ^ J d'^kFk{t)ala^k + F^{t)aka-k 

+Ekit) {2alak + S^''\0)) , (90) 
where the factor Ek {t) for the diagonal terms is, 

Ekit) = \vk{t)\'' + Wlit) + l/A[Ukit)/Ukit)]''}\vk\'' 
+ l/2Uk{t)/Ukit)dt\vkit)W (91) 

and the off-diagonal term Fk{t) is given by 
Fk{t) = {vk{t)f + {iul{t) + l/i[Uk{t)/Uk{m{vkf 



l/2Uk{t)/Uk{t)dt{vk{t)f 



(92) 



Therefore there are no time-independent eigenstates that rep- 
resent the physical vacuum state, one that minimizes the ex- 
pectation value for the Hamiltonian, 

((t)0|7?fc(t)|(,)0) = J j d-'kEkit), (93) 

for all times. 

Yet somehow, in order to obtain a meaningful statement 
about the particle production in our effective curved space- 
time, we need to deal with or circumvent this problem. 

a. Instantaneous Hamiltonian diagonalization: One 
possibility to find an approximate vacuum at an instant of 
time — say t — — is to define the vacuum state |(to)0) of 

the instantaneous Hamiltonian H{tQ), where we project our 
mode functions instantaneously onto a plane-wave basis with 
the dispersion relation ujk{ta)- Within this approximation we 
find for the coefficients in Eq. &U[: 



(94) 
(95) 



Ek = \vkiW + i^Uto)\vkr; 

Fk = {vk{t))^ + ^l{to) {vkY ■ 



It can be shown that the expectation value for the instanta- 
neous Hamiltonian — up to an arbitrary phase — at t = to 
is minimal for mode functions Vk — \ / uJk{to) and Vk — 
i ^/uj ^ iujk Vkito). 

Note that at to the two coefficients Fk{to) = and 
Ek{to) = 2ujk{ta), and that under this condition the instan- 
taneous Hamiltonian (see Eq. (90 1) is diagonal. Therefore it 
is referred to as the vacuum state of instantaneous Hamilto- 
nian diagonalization. 
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In a time-dependent problem the vacuum state changes, and 
therefore one has to give up on this particular definition of vac- 
uum for times t > t^. The association between the \(3k\'^ Bo- 
goliubov coefficients (as outlined above) and physical quasi- 
particles requires that we approach — either asymptotically or 
abruptly — flat initial and final state spacetimes. Only a time- 
independent spacetime perpetuates its vacuum state, where 
LOk{t) = LOk, and thus the normalized mode functions 



1 



Vk {t) = , — — exp(i Wfc t) , 
\/wfe(t) 



vl{t) 



exjp{-iujk t) , 



(96) 
(97) 



in this case represent a physically meaningful vacuum state 
for all times t. 



Bogoliubov transformation in emergent spacetimes 



The Bogoliubov transformation outlined in Sec. Ill A 3 can 
be applied to time-dependent problems, as long as we avoid 
making any statement about the presence of "real" quasi- 
particles unless the expansion has started from an initial time- 
independent spacetime, and eventually approaches a time- 
independent spacetime. We also require uj\ > 0, such that at 
the beginning and the end the mode functions are simple (nor- 
malized) in- and out-going plane waves. For a constant, but 
imaginary frequency, ujf < 0, the instantaneous Hamiltonian 
is still diagonal, but its expectation value no longer has a min- 
imum. (A detailed treatment of this problem in spacetimes 
emerging from a Bose gas can be found in Ref. |9|. There 
the underlying physical model exhibits a negative scattering 
length.) 

To apply the Bogoliubov transformation (see. Sec. |III A 3| 
and see Eqs. ( [86p7| ) to our specific problem it is necessary 
to consider that the original problem was defined in terms of 
two linearly independent field variables, the field operator 6, 
and its conjugate momentum 11^. The connection conditions 
arise from the necessity that the field operator, 

= lim|^(t-e)-0(t + e)| = 0, (98) 

and its conjugate momentum (on the emergent spacetime). 



= lim 



i{n,^(i-e)-n^(< + e)} = 0. 



(99) 



have to be continuous at all times. (A more detailed treatment 
of this problem in the context of emergent spacetimes 
from Bose gases can be found in |9|.) Note, that these 
two conditions have to be fulfilled for any arbitrary (albeit 
physically reasonable) change in the contact potential U{t). 



In Sec. Ill A 1 we have already pointed out that our emergent 



spacetime is special in the sense that all time-dependence 
(in our model) has to be implemented via U{t), which also 
shows up in the overall conformal factor in the equation of 



motion, see Eq. (66i. 



a. Bogoliubov coefficients for a time-dependent emergent 
spacetime: Let us consider a contact potential defined as 
follows, 

Uk{t) = gk{t) easit - to) + hk{t) Q^isih - t) . (100) 
We are dealing with different sets of mode functions in each 



region; 



{vk,vl) for t < to; 



and — 



{uk,ul) for 



t > to- (Here 0hs is the symbol for the Heaviside step func- 
tion.) For the time being we do not make any further assump- 
tion except that the mode functions have to be a solution of 
the equation of motion ( [66] l in their respective region. In par- 
ticular, we do not yet assume the modes are normalized. The 
connection conditions, given in Eqs. (98 i and (99i, in com- 



bination with the Bogoliubov transformation — see transfor- 



mation matrix (83 1, mapping between two pairs of complex- 
conjugate mode functions (it is always possible to choose such 
solutions) — provide us with two matrix equations (that is, 
four component equations). 




Mdt 




(101) 



, (102) 



for the four unknown transmission coefficients ak, a^., (3k, 
and PI — which are contained in the matrix M. Fortunately, 
due to our specific choice of complex mode functions, we only 
need to calculate ak and /3fc, and get the others by calculating 
their complex conjugate. We obtain 



ak{t) 
and 

Pk{t) 



2 {g UkV*k - h UkV*k) + vluk [gj^-^h 
2^^W[uk,ul] 



(103) 



2^ W[uk 



(104) 



A somewhat time-consuming, but trivial calculation shows 
that further conditions are necessary to obtain a Bogoliubov 
transformation in the desired normalized form (|84li, since in 



general ( 103 1 and ( 104 1 lead to 



\ak{t)\'^\l3k{t)\'^ 



W[vk,vl] 
W[uk,ul 



y^l 



(105) 



Hence it is necessary to choose a consistent normalization 
condition, iy[wfc,wy — W[uk,ul.] = constant, so that 
|afe(i)p - |/3fe(t)P = 1. Also note that for multiple- 
step events, meaning multiple Bogoliubov transformations de- 
scribed by the compound matrix Af„ • il/„_i • • • M2 ■ Mi, the 
physics is independent of the particular choice for the normal- 
ization of the intermediate mode functions, since 



det(M„ • • • Ml) = det(M„) • • • det(Mi) = 



(106) 
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Our main focus in this article is to calculate the particle 
production in a FRW-rainbow metric, and therefore we wish 
to restrict our problems to cases where we start and end in a 
physically meaningful vacuum state. This way we are able 
to connect P with the occupation number for real conden- 
sate excitations in an effectively expanding spacetime. In our 
companion paper, see (T\, we investigated various expansion 
scenarios. We will now concentrate on the most relevant ex- 
pansion scenario in terms of cosmology, i.e., de Sitter like in- 
flation. 

There are three principal different cases, in terms of U{t), 
that are of interest when addressing the problem of inflation 
in emergent spacetimes: 

(i) U(t) discontinuously, connects two flat spacetime regions: 

The sudden case is indirectly interesting for inflation, 
since the de Sitter expansion approaches the sudden 
case for infinitely fast expansion |34|. But, we will 
show shortly that in emergent spacetimes — due to the 
nonperturbative ultraviolet corrections — the particle 
production for the extreme limit is finite, thus the Bo- 
goliubov transformation for the de Sitter case is always 
well-defined. We discuss this in Sec. 1111 B I 

(ii) U(t) continuously, but not continuously differentiably, 

connects two flat spacetime regions with a flnite de 
Sitter like phase in between: This calculation can be 
carried out in the hydrodynamic limit, but fails to be 
a good approximation for infinitely long-lasting infla- 
tion in emergent spacetimes. Due to an unsolvable sec- 
ond order differential equation in the eikonal limit, we 
present a qualitative analysis, that should be compared 
with the numerics presented in our companion paper 
n. We discuss this in Sec. |lVCl| 

(iii) U{t) is a smooth function everywhere, such that we are 

left with one de Sitter like region: The qualitative analy- 
sis shows that in our particular emergent spacetime the 
de Sitter expansion has, in the infinite past and in the in- 
finite future, two distinct physical vacua, and therefore 
we are able to predict the existence of time-independent 
real (unambiguous) quasi-particles created during an in- 
finitely long-lasting expansion. This might at first seem 
of less interest for the condensed matter community, but 
as our simulations show, this can be realized for a suf- 
ficiently long expansion time. We will see that for the 
right parameter choice the particle spectrum approaches 
a characteristic final shape, for which we can numeri- 
cally determine the form of the final particle spectrum. 
We discuss this in Sec. llVC4] 



The case of a particle production for a sudden transition has 
been previously explored by Jacobson for a parametric oscil- 
lator ||35]| . (Though the underlying physics is rather different, 
there is also a model for sonoluminescence that is based on a 
rapid change in refractive index — that model shares many of 
the mathematical features encountered in the present calcula- 
tion |l2l]|36l|37l[38l[39l|40l|4T].)A more closely related work 
on sudden changes — between Lorentzian/ Euclidean signa- 
tures — in emergent spacetime has recently been carried out 
int9J. 

To calculate the particle production in the limiting case 
of a sudden expansion, we consider the situation where the 
atom-atom interaction is "instantaneously" switched from U 



to at some time Iq. The scale function (58 i is given by 



m 



1 - ^ ) Qnsit - to) . 



(107) 



This corresponds to a change between two regions with dis- 
tinct dispersion relations, connected at to. We choose as posi- 
tive and negative frequency mode functions for t < to. 



Vk{t) = — 1=- and vl{t) = — -=■ 



(108) 



and for t > to, 



Uk{t) = 



t 



/, ,ou 



and Mfe(t) 



(109) 



these are normalized using Eq. (|82]). As explained in 



Sec. Ill A4a these mode functions only represent physically 
meaningful vacua in flat spacetime regions. 

Straightforwardly, we can apply Eqs. ( 103 i and ( |104| i to 
calculate a and /3. We get for the Bogoliubov coefficients; 



ak{to) 
h{to) 




(110) 



('^r'+"") (111) 



Using the transformation laws between the in and out Pock 



state operators we obtain the occupation number density (88 1 
of out -particles in the in-vacuum. 



„out 



1 



(112) 



B. Finiteness of particle production in emergent spacetimes 

In a universe that is subjected to a (finite-size) expansion, 
and in particular a sudden variation in the size of the uni- 
verse [34 1, there is a relatively simple way to calculate an 
upper bound on the particle production. 



where we assumed that the eigenfrequencies are real. Note 
that this relation is valid beyond the eikonal limit, such that 



we are dealing with a nonlinear dispersion, see Eq. ( 59 1, in 
both regions; 



out 



4Vi + ^V^', 
4Vi/^ + fcV^' 



(113) 
(114) 
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Thus we get 



(115) 



It is perhaps surprising how simple the final result is, despite 
the additional technical machinery required to derive it. Of 
course, particle production is momentum-independent only 
within the hydrodynamic limit, when 



in/out 



in/out 



(116) 



The number density in the infrared limit, when k < K, can be 
simplified to 



1 



^hydio 



^1/4 „ ^-1/4 



(117) 



and is in agreement with the results obtained in "ordinary" 
momentum independent spacetimes with sudden jumps in the 
dispersion relation, instead of the scale factor; for example 
see |35|. Bogoliubov transformations with infinite total par- 
ticle production per unit volume are ill-defined, as the out- 
vacuum state loutO) can be written in terms of a normalizable 
linear-combination of in-particle states, only if |/3fep 
faster k^'^ for large k. To see an explicit derivation consult, 
for example, Ref . B2l . We shall now show how due to non- 
perturbative corrections in our FRW-rainbow spacetimes, this 
problem does not appear Here the total number of particles 
produced in the ultraviolet limit remains finite under any cir- 
cumstances. At large momenta 



lim Uk 



X2 



K 
2k 



0(1/^*^), 



(118) 



and separately we see that for an infinitely large expansion 
that 



1 + fc2 _ fc)2 

lim fik — , : — 

x^oo 4 k y/K^ + fc2 



(119) 



We can also combine these two limits, and obtain Uk 
[if/(2 /c)]"*. The quantum pressure term supresses particle 
production at high momenta — effectively because at high 
momenta the quasi-particles are free and do not "see" the 
evolving spacetime. (The high-momentum limit of the rain- 
bow metric is "static".) 

As a consequence the number particle production per 
spacetime volume. 



N. 



<k 



dk k" 



r d-4 



(120) 



for the here relevant cases of 2 or 3 spatial dimensions, is fi- 
nite. 



It is also possible, using Eq. ( 119 1, to exactly calculate the 
total number of particles produced during, cq cq/VX, for 
X ^ oo, an infinite expansion: 



Nr. 



(121) 



which is also finite. 

Particle production in our emergent spacetime is a good ex- 
ample as to how corrections — at the level of effective field 
theories — from the underlying microscopic structure can cir- 
cumvent the problem of an infinite number density. On the 
other hand, it seems naive to assume the validity of common 
quantum field theory on all scales. If spacetime is indeed the 
infrared limit of a more fundamental theory, the presence of 
ultraviolet deviations seems to be unavoidable. 



IV. FRW-RAINBOW SPACETIMES IN (2 + 1) 
DIMENSIONS: 
CONTROLLED PARAMETRIC EXCITATIONS IN 
2-DIMENSIONAL CONDENSATE 

In the following we investigate the robustness of the cosmo- 
logical particle production process against model-dependent 
ultraviolet corrections in the dispersion relation. The modifi- 
cations arise from the microscopic substructure, therefore are 
of non-perturbative nature and not the result of some pertur- 
bative loop-calculations. Related work on Hawking radiation 
from acoustic black holes has been carried out, where the var- 
ious authors mainly studied the problem in media with sub- 
sonic and supersonic dispersion relations lfT5l[T6l . In addition, 
similar efforts have been made in the field of "conventional" 
cosmology, by implementing ad hoc trans-Planckian modifi- 
cations to the system, for example see Il43l l44l l45l . This mo- 
tivated us to extend our main intention — testing the analogy 
— to a detailed study of the deviations (from standard cos- 
mology) in our specific model. We demonstrate this with a 
particular physically relevant model, a rapidly exponentially 
growing effective universe, how "conventional" and "mod- 
ified" models naturally merge in the infrared limit at early 
times. Whereas, they show significant deviations at late times, 
yet for a sufficiently short-time expansion, the resulting parti- 
cle spectra — in systems with linear and nonlinear dispersion 
relations — are very similar to each other. 



A. (2 + 1) dimensional FRW-rainbow spacetimes: 
2 dimensional condensates with time-dependent atomic 
interactions 



In Eq. ( 62 1 we established the concept of a FRW-type rain- 
bow universe to describe the (2 + 1) or (3 + 1) dimensional 
spacetime (that is, for d = 2 and d = 3 space dimensions 
and 1 time dimension) as experienced by small quantum fluc- 
tuations in a 2 or 3 dimensional condensate. As pointed out 
in Sec. II B besides the dimensionality, the scale factor h{t) 
is also different for 2 and 3 dimensions. For example, if we 
wish to mimic a de Sitter geometry, where the scale factor in 
the FRW Une element, see Eq. ([TJ, takes the form 



a(T) = cxp(iJ t) 



(122) 



we have to choose the scale factor for the atomic interactions. 



exp(—t/ts) 

ts/{t + ts) 



for 
for 



(123) 
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A more detailed treatment of this problem — that originates 
in the discrepancy of laboratory t and proper time r for the 
3-dimensional case — can be found in [3|. In 2 spatial dimen- 
sions — where a{2) = — the effective line-element (62 1 is 
then given by 



dsL2 = 



(124) 



here laboratory time t and proper time r are of the same form, 
and we obtain the very simple relationship, 



ak{t) = bk{t) 



-1/2 ^ 



1 



(125) 



between the two scale factors. To connect tg, the scale unit in 
the laboratory, to the Hubble parameter H, we need to choose 



H 



1 

2U 



(126) 



We will soon revisit the Hubble parameter, when we discuss 
its physical motivation, and how the Hubble parameter applies 
to our de Sitter-rainbow spacetime. For practical reasons we 
shall be mainly interested in (2 + 1) spacetime dimensions, 
simply because the numerical simulations we wish to com- 
pare the theory with are more easily carried out in 2 space 
dimensions. 

We will always use a{t) to refer to the FRW scale factor; if 
we ever really need the scattering length we will refer to it as 

I3.scatt(^)- 

In what physical respect do we really have an "expanding 
universe", given that the condensate is physically contained in 
a fixed volume V7 A decrease in the scattering length corre- 
sponds to a decrease in the speed of sound propagating in the 
condensate; therefore any quasi-particle excitation will prop- 
agate with decreasing speed in the condensate as time passes. 
To an "internal" observer at rest in the effective spacetime, 
and communicating by means of acoustic signals, a decrease 
of the speed of sound is indistinguishable from an isotropic 
expansion of the spatial dimensions. In contrast "external" 
observer made out of real particles belongs to the microscopic 
world, and thus is inevitably bound to laboratory frame and 
not part of the analogy | 



B. Inflation in emerging universe: 
Exponentially decreasing atomic interactions 

Maybe the most interesting cosmological case to study in 
our emergent spacetime is the de Sitter universe, where the 
scale factor is given by an exponentially expanding (or con- 
tracting) universe. The concept of inflation was introduced si- 
multaneously around 1981 and 1982 by Guth |47 1, Linde |48 1, 
and Albrecht and Steinhardt |49| to explain the homogeneity 
of the temperature observed in our universe, beyond casually 
disconnected areas. Not long after {e.g., see Guth, Hawk- 
ing |50|, Bardeen ISTI . Turner 15211 and Brandenburger |53|) it 
has been realized that inflation also accounts for the existence 



of the perturbations in our universe today. These perturbations 
in the form of a slight deviation from a uniform temperature 
in our cosmological microwave background (CMB) has been 
measured in 2001 by Netterfield et al. To date, inflation — a 
phase of a rapid de Sitter-like expansion between a predating 
radiation a{t) ~ t^/"^ and postdating matter a{t) ^ t^^^ dom- 
inated area — seems to be the most plausible explanation for 
the CMB map. For more details see |54 55 1. While we expect 
deviations from the standard picture in our emergent rainbow 
spacetime, it seems desirable to hold on to the concept of a 
thermal spectrum resulting from an exponentially expanding 
universe. We will come back to this point at the end of this 
section and show that this hope might not be in vain. 

Before we focus on the particle spectrum in our de Sitter- 
rainbow metric, we would like to introduce some physically 
important parameters. They will be of great value throughout 
the remaining part of this article, and are necessary to under- 
stand the particle production process for a 2d analogue FRW- 
type universe. Let us start with the obvious question: 



1. Why rainbow spacetimes? 

The modified rainbow spacetimes owe their names to their 
momentum-dependence. This kind of modification can be 
absorbed into the time and (now also) momentum-dependent 
scale factor For a de Sitter like universe (see Eqs \\22\ and 
( 123 1) in the hydrodynamic limit, we get 



ak{t) 



^exp{-2Ht) + {k/KY 



(127) 



for the modified scale factor ( |122[ l; in 2 spatial dimensions. 
Thus the hydrodynamic. 



crossover. 



and free particle. 



exp{-2Ht) > |fc/if|2 



exp{~2Ht) - \k/K\^ , 



exp{-2Ht) < |fc/if|2 



(128) 



(129) 



(130) 



limits are a matter of dividing the spectrum into appropriate 
energy regimes at a particular time t. It is interesting that for 
early times — when the interactions between the atoms are 
strong — we naturally approach the hydrodynamic case. 



lim ak{t) a{t) , 



(131) 



in the sense that most modes are phononic, and therefore 
larger and larger fc-values are covered by "conventional" 
FRW-type quantum-field-theory. 

Quite the contrary occurs after an infinitely long-lasting ex- 
pansion, where all modes behave as free particles. 



lim afc(t) 

t — *+C30 



\K/k\ , 



(132) 
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and the universe, as seen by a mode with the wavelength k, 
will approach a final finite fixed size. 

Due to this fundamental difference between our analogue 
model and the "theory" we wish to mimic, we know already 
that there will only be a finite time-period — its length de- 
pends on K, and therefore on the tunable initial interaction 
strength U{0) — Uq — beyond which the analogy breaks 
down. Note that the particle production process must natu- 
rally come to an end, when the expansion rate slows down to 
zero. That the effective expansion rate of the FRW-rainbow 
universe approaches zero will be shown next, but before we 
would like to answer the question posed in the headline of this 
subsection, with a simple illustration for the de Sitter-rainbow 
scale factor in Fig. [T] There we plot the emergent rainbow 
scale factor a^it) for each k mode using different colors — 
gradually changing from dark red for infrared modes to dark 
blue for ultraviolet modes. The resulting color-spectrum is 
reminiscent of on the color spectrum obtained from real rain- 
bows. 



2. What about the Hubble parameter? 

The de Sitter universe is special in the sense that its Hubble 
parameter, the "rate of expansion" 



H 



d(t) 
a{t) 



is constant. The de Sitter universe is a solution of the Einstein 
equations with a positive cosmological constant. A; i.e., H oc 
a/a. If the acceleration of our universe can be put down to 
this cosmological constant, the universe will expand forever 
— and further dilute the matter and radiation distribution in 
our universe — until it approaches the de Sitter spacetime ll56l 
|57J. 

What is the situation in our emergent de Sitter universe? 
The rate of size change in the emergent de Sitter universe us- 
ing [127] is given by, 

H, = H exp(-2g^) 
* exp(-2i/i) + (fc/if)2 ' ^ ' 

a momentum-dependent rainbow Hubble parameter. At early 
times, or for phononic modes, when ak{t) — > a{t), the rain- 
bow Hubble parameter, 

lim Hk{t)^H, (135) 



reduces to the conventional Hubble parameter, while for late 
times, 

HK'' (exp(-fft)/fc)2 



lim Hk{t) 

t — >-+oo 



0. 



(136) 

stops expanding 



The universe gradually — mode by mode 
as modes leave the phononic regime. 

Therefore, in our particular scenario the long-time- 
kinematics of the particle production process is determined 
by the non-perturbative corrections from the substructure, and 
not by the emergent spacetime picture. Any speculations with 
respect to an everlasting expanding universe are, (within our 
model), in vain. 



3. Characteristic value for quantum process? 

To get a grasp on the cosmological particle production 
process, these two parameters (the effective scale function 
at {t), and the effective Hubble parameter Hk (t)) are not quite 
enough. They are a good measure to describe the kinemat- 
ics of the emergent gravitational field, but one also needs to 
know how the microscopic corrections affect the energy of 
the modes. Only the ratio between the mode frequency. 



LUk{t) = ujo Vexp(-2iJt) + {k/Ky 



(137) 



and the Hubble frequency can tell us whether the mode will 
be disturbed by the classical background or not. This result 
has been established for quantum field theory in conventional 
de Sitter spacetimes, and we review this point in the following 
section. 

For now, we simply transfer the qualitative description 
known from conventional cosmological particle production, 
to our rainbow spacetimes. In this spirit we define the ra- 
tio between the modified dispersion relation and the effective 
Hubble parameter as follows: 



7^fe(<) = 



ujk{t) _LUo {cxp{-2Ht) + {k/K)^f/^ 



exp{-2Ht) 



(133) Within the hydrodynamic limit this ratio simplifies to. 



Rk{t) 



H 



exp{—Ht) 



(138) 



(139) 



a monotonically decreasing function with time. Again, we 
calculate limits for very early. 



lim 7^fc(^) 

t — * — oo 

and very late times. 



lim 7^fc(^) 

t-^-\-OQ 



H 



UJQ 

H 



cxp{-Ht) = Rk{t), (140) 



{k/Kf exp{+2Ht) , (141) 



to see once again, that for late times /trans-phononic modes 
the deviations in our emergent spacetime play an important 
rule. 

We will show that this ratio is sufficient to understand, and 
therefore predict qualitatively, the particle production process 
in our analogue spacetimes. 



C. Quantum field theory and rainbow inflation: 
Excitations from an exponentially clianging interaction strengtii 



The equation of motion in the presence of time-dependent 
atomic interactions is given in Eq. (|66ll. In the de Sitter rain- 



bow universe, with the scale factor ( 127 1, and the dispersion 



relation (137i, the harmonic oscillator frequency (70 1 simpli- 
fies to 



nut)' 



' bk{t) - H' + Ak{t)' 



(142) 
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(c) Scale factor quantum pressure effects ; is = 1 X 10 ^. (d) Scale factor quantum pressure effects; ts = 1 X 10 ^. 



FIG. 1: (Colors online only.) In this figure we plot the logarithm of the scale function ak{t) for each fc-value — for k G [9, 191] — in 
a different color. The different colors encode the energy of the modes: Gradually changing from low-energy / infrared (dark red) to high- 
energy / ultraviolet (dark blue). The upper row shows the behavior of the scale function in the hydrodynamic limit. While the rainbow-scale 
function — shown in the lower row — approaches the hydrodynamic limit for low-energy modes, the ultraviolet modes show strong deviations. 
Note, that in the infinite past all modes are phononic, and therefore ak{t) a{t). The black dots indicate the time-dependent crossover 
(phononic to trans-phononic) in every quantum mode. Parameters are CNhit — Q) = 2 x 10^, A'^o ~ 10^ and X = Ax 10®. 



where the effective (or rainbow) scale factor is given by 

hk{t) = exp[-2i?i] + {k/Kf , (143) 
and the last term is given as 



. .,^2 ^ rr2 [4exp(-2H^) + (fc/X)^] {k/Kf 
' [exp{^2Ht) + {k/KYY 



(144) 



which can be neglected in hydrodynamic approximation. 
Therefore the equation of motion for our rainbow de Sitter 
spacetime is given by a rather complicated differential equa- 
tion, where common techniques for solving second-order dif- 
ferential equations fail. 



However, this is not the end of the story, since at early 
times, or for phononic modes we get. 



bk{t) ^ b{t) 



lim 

t — > — oo 



where 







(145) 



b{t) = exp(-2i7<) , 



(146) 



and within this limit the equation of motion simplifies to a 
differential equation solvable by Bessel functions. Within this 
regime, our emergent spacetimes can be regarded as a good 
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toy model for quantum field theory in conventional de Sit- 
ter spacetimes. However, we need to apply special boundary 
conditions, that is to set limits on the time-duration and the 
fc-range. Only then will the cosmological particle production 
process be robust against the model-specific non-perturbative 
corrections that are present in a spacetime emerging from a 
gas of Bosons. 

As the expansion continues, the more modes will cross-over 
from phononic to trans -phononic, and finally for an infinitely 
long-lasting expansion, we get 



lim 

t — ^+00 



bk{t) 



{k/Kf 



H 



{k/Kf 



(147) 

and therefore are left with no phononic regime whatsoever 
That is, in the infinite future all excitations of the fluctuations 
of the system are free-particle like. In this regime the quantum 
fluctuations are decoupled from the collective behavior, the 
emergent spacetime picture, and thus are static. 

Therefore, without any further considerations, we de- 
duce that in the infinite future we will always end up with 
well-defined out-states, that represent the vacuum state of 



instantaneous Hamiltonian diagonalization, see Sec. Ill A 4 a 



The normalized mode functions in the infinite future are given 



by Eq. (lK)9|l, where =U()\k\/K. 



We will subsequently revisit this point in Sec. IV C 2 but 
for now our concern will be to focus on a regime will be to fo- 
cus on a regime — the hydrodynamic limit — where the prob- 
lem of particle production is tractable, and there is a straight- 
forward description of the physics. 



1. Toy -model for conventional inflation 

To apply the method of instantaneous Hamiltonian diago- 
nalization to particle production from a finite de Sitter phase, 
we assume a continuous function for all t, but at the times 
and tf the function is not differentiable: In the hydrodynamic 
limit 



hk{t) b{t) 



(148) 



we assume 



b{t) = l+[b^^{t)-l] eHs(i-io)-[fo''(<)-iA] eHs(i-t/); 

(149) 

compare with Fig. Q. 

Here h'^^{t) = exp[-2i/(t - io)], such fliat 6(to) = 1, 
and b{tf) = 1/X. In addition, we define LUk{to) = ^'k^ ™d 
Wfc(i/) = 



such that for t < to we choose the in-mode 
functions, given in Eq. (108 1, and similarly for t > tf we 



choose the out-mode functions, given in Eq. ( |109| l, for the 
initial and final flat regions. 

Note that for any laboratory set-up the time-scales are 
finite, but given that this calculation is only valid for modes 
that are of phononic nature both before and after the expan- 
sion, we are not free to choose the expansion time arbitrarily. 




FIG. 2: Time-dependence for atom-atom interactions to mimic a 
de Sitter universe in a 2 dimensional Bose-Einstein condensate. 



The Bogoliubov transformation is slightly more difficult 
here, where two single-step processes have to be carried out: 

Vkit) = M{tf) ■ vf{t) = M(tf) ■ Mito) ■ Mt) ■ (150) 



For the mode functions during the de Sitter phase, which 



are v'^^ and {v/. 



dS 



)* for t G [io,t/], we solve the harmonic 
oscillator equation with the time-dependent frequency given 



in Eq. ( 147 1. The solution is a linear combination of first or- 
der Hankel functions of the first H^^' and second H^^'' kind. 
The normalized mode functions — using the Wronskian con- 
are given by 



dition (82 1 



(151) 
(152) 



Note that in this representation the mode functions are indeed 
a set of complex conjugate functions, since (H^^^)* — H^^^^; 
see for example |58|. The argument of the mode functions 
is Rk{t), the ratio between the mode frequency and the Hub- 
ble frequency, see Eq. ( |139| l. These mode functions are only 
valid within the hydrodynamic limit, i.e. for modes that satisfy 
\k/K\<$:exp{-Ht). 

To calculate the Bogoliubov coefficients at each step we 
again apply Eqs. ( 103 1 and ( 104 1. Below we will briefly out- 



line the necessary steps to calculate the final Bogoliubov co- 
efficients for the whole process. To begin with we write down 
the single-step 

Uit) = git) emit - to) + hit) e^ito - t) . (153) 
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(a) t = 0.25 X tf;ts = 5 X W'^. 
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(c) t = 0.75 xtf;ts = 5x W'^. 
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(b) t = 0.50 xtf;ts=5x 10'^ . 




100 150 200 

Ikl 

(d) t = 1.00 xtf;ts=5x 10-5. 



FIG. 3: The four figures show time-slices of the "quasi-particle" production Nk in a de Sitter like spacetime, for ts = 5 x 10~^. The blue dots 
represent the actual data obtained from the simulations for a realistic Bose gas with time-dependent interaction-strength |T |. The dashed green 
line shows the theoretical results for a finite de Sitter calculation obtained in the hydrodynamic limit, as presented in Sec. |lV C"T| The vertical 
(red) line indicates the crossover between the two relevant regimes, the excitations being phononic on the left and trans-phononic / free-particle 
like on the right. 



Bogoliubov coefficients, 

ak{t) = ^ (2W[uk,vl] + vluk{h/h- g/g)^ (154) 
Pk{t) = ^^(2W[vlul]-vluUh/h-g/g)), (155) 

appropriate to a continuous, but not continuously differen- 
tiable scale function at the time t. At the first step, where 
we map from plane waves onto Bessel functions, and for our 
specific scale factor before the step, g{t) — 1, and after the 



step, h{t) = exp[— 2ff(t — to)], we get 



ak{t) = -{W[vk,{vfr]~Hvlvf} , (156) 
Pk{t) = l-^{W[vl,{vfr] + Hvl{vfr} . (157) 

The final coefficients for our specific mode functions are given 
by 



(158) 
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and 



7^ 

i2 



{H(2'(i?l")-zHf)(i?l")} . 

(159) 

In the last two equations we have used R]^^ = Rk{to) and 
cj" = ujk{to)- Instead of calculating the Bogoliubov coeffi- 
cients for the second step explicitly, we suggest a little short- 
cut. The transformation at the first step was a mapping from 
plane waves, represented as Vk and w^, onto the de Sitter mode 
functions v'^^ and (u^^)*: Vk = M{t)vf^. At the second step 
the calculation is exactly the opposite, a mapping from the de 
Sitter mode functions onto plane waves: v^^ = M{t)^^Vk- 
The inverse of the transition matrix displayed in Eq. ([83]l is 



(160) 



where we have employed \ak\'^ — |/3fcP = 1, since we map 
between normalized mode functions. Altogether we can for- 
mally derive the Bogoliubov coefficients at the second step 
as; 



> [a 

in — ^out 



(161) 



pntf) ■■ /?i"(io)^^-[/3r(^/)], (162) 

in — >out 

Thus we can write down the Bogoliubov coefficients at the 
second step, at i = tf, without any farther calculation being 
required, as 



and 



2i 



k g"^'^ 



(163) 



/_ DOUt , - 



Herei?^"* = i?fe(t/) and( 



i^kitf). 



(164) 



The overall particle production, at later times t > tf can be 
obtained by simple matrix multiplication, since 



a 



final 



final 



Mt>t, = M{to) M{tf) = y ^•^t.xy (^finai)* j -(165) 
and therefore we get for the final BogoHubov coefficients 



final 



„ in „ out I Qin { /OOut\ * 

"fc "fc + Pk [Pk j 



(166) 



/□final „ in pout i pin ( „ out \ * / 1 iin\ 

Pk = f^kPk +Pk ("fc ) • (167) 

Note that these coefficients can in some sense be considered 
as time-dependent Bogoliubov coefficients, where 



final(t') 



+p'k{pT^'y, (168) 



^finaKO ^ ^n.^out(t')+^in(^out(t'))*^ (Igg) 

where we project at any instant of time onto a plane wave 
basis at a particular t' with the eigen-frequency ujk{t')- But, 



as we pointed out in Sec. Ill A 4 as long as the expansion 
is continuing, the corresponding out -frequency modes do 
not represent a physically meaningful vacuum state. Nev- 
ertheless, for t > tf, we stop the expansion and force the 
Hamiltonian of the system to become static, so that we are 
able to associate n^'^^^ — |/3f"^'p with the "real" mode 
occupation of a mode k per unit volume. 

A lengthy but straightforward calculation now yields the 
mode occupation number after stopping the de Sitter like ex- 
panding phase. The number of particles produced depends 
only on the initial and final values of the frequency ratio 
Rk{t), respectively and as: 



64 



oin DOut 

■Kfc Rk 



/ Din^ tt(2) / poutN 

[HW(4")H(^)(i?r) 



.out) _ jj(2)(-^i. 



jj(l)(^out)]2 



H(^)(i?i")HW(i?r)]^ 



joutM2 



H«(i?r) 



(2), 



joutM2 



2Hf)(i?r)H^'^(i?nx 
[H«(i?i-)HW(i?r) 



H«(i?L")H«(i?r) 



2Hf)(i?I")H^^^(i?r) 

[H«(i?i")HW(i?r) 



(2), 



H 



(^Hi?r)H«(4")] 



(170) 



At first sight this formula seems to be rather complicated, 
but there are two relatively simple and feasible consistency 
checks. 

For example, we expect no particle production in the limit 
where the universe has not changed at all. It can easily been 
seen that for 



lim n 



dS 
k 



0, 



(171) 



the occupation number indeed goes to zero. 

Also interesting are limits resulting form an infinitely slow 
(adiabatic), that is iJ ^ 0, and an infinitely fast (sudden) 
H ^ oo expansion. In either case we replace the ratios with 



: u'^/H and 



7i7, hold and fixed. 



while we allow the Hubble parameter H to vary. Such a pa- 
rameter choice requires Ht^xp — constant, and therefore we 



expect the expansion time t^ 



= t 



f 



to to behave inverse 



proportional to H. The leading order, obtained from a Taylor- 
series-expansion around H = oo, 



H 



lim n'jf 



,out ) 2 



0{l/H^), (172) 



is as expected — since (tf — to) — in agreement with the 
result from the sudden calculation within the hydrodynamic 
limit, see Eq. (112i. In contrast, a Taylor-series expansion 
around H — yields 



lim nf^ 



0. 



(173) 
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As expected, within the Umit of infinitely slow expansion, we 
produce zero particles, and we recover adiabatic invariance. 
Furthermore if we consider the asymptotic expansion where 
1 <C -RJf ^ V^, and employ the asymptotic limits of the 
Hankel functions [58.1 . we get at linear order: 



nk 



H 

2tt 



(174) 



Of course, we cannot rely on our calculation in the case of 
an infinitely long-lasting expansion, since it is based on the 
validity of the hydrodynamic limit, which is completely inap- 
propriate for t ^ ts- However, in our previous numerical sim- 
ulation of a "realistic" Bose gas 1 1 1 we compared our theoreti- 
cally obtained result with short-time expansion scenarios, and 
were able to match them to the phononic part of the particle 
production spectrum, see Figs. [3] (In [T] we applied the trun- 
cated Wigner method to a Bose-Einstein condensate, where 
we selected a nonlinearity and atom number such that most of 
the modes were phononic at the start of the simulation, and 
that for our choice of parameters the effects of back-reaction 
and mode-mixing were minimal.) 

In Figure [3] for = 5 x 10~^ we see an excellent agree- 
ment between our theoretical predictions (dashed green line), 
and the numerical data (blue dots) for the "quasi-particle" pro- 
duction in the phononic regime (left of the vertical red line). 
As shown in Fig. [3] at the end of the expansion we see that 
almost all excitations are trans-phononic, hence the analogy 
eventually breaks down. But, to answer the question stated at 
the beginning of this section, where we asked about the appli- 
cability of our specific model to mimic cosmological particle 
production: Yes, the BEC can be used as an analogue model 
for cosmological particle production within certain limits, as 
our numerical results clearly confirm. 

The calculations we have presented are somewhat tedious, 
but they are more than worth the effort since there is a fun- 
damental lesson to be leamt from Eq. ( 170 1: It establishes 



our previous intuition (see Sec. IV B 3 1, that the characteristic 
value controlling the particle production is the ratio, Rk (i) be- 
tween the mode frequency, and the Hubble frequency. Com- 
parison with Eq. ( |170| l shows explicitly that the final quan- 
tity of particle production — for truly phononic modes be- 
fore, during, and after the expansion — only depends on the 
initial and final i?™' frequency ratio. This motivated us 
to extend the role of the frequency ratio beyond the hydrody- 
namic limit, and compare our predictions with the data ob- 
tained from our simulations, where the non-perturbative cor- 
rections are included; see yj for details of the numerics. 



2. Qualitative behavior of quantum fluctuations 

The process of cosmological particle production in an ex- 
panding/collapsing universe can be qualitatively understood 
in terms of a single parameter, the frequency ratio TZk{t) as 
given in Eq. 



138 1 



First, we explain the connection between 
the qualitative behavior of the particle production process and 
this frequency ratio ( 139| l in the hydrodynamic limit, then we 
transfer these ideas to the emergent rainbow metrics we have 



introduced in Sec. lIIBI 

a. Qualitative behavior of particle production in the hy- 
drodynamic limit. There is a relatively simple way to under- 
stand the qualitative evolution of mode functions in an expo- 
nentially changing universe, by considering out the effective 
harmonic oscillator equation for the auxiliary field. Within the 
hydrodynamic limit we get 



(175) 



where we used the equation of motion, see Eq. ( 66 1, for the 



effective time-dependent harmonic oscillator frequency for 



phononic modes given in Eq. ( 145 1. 

Above we have shown that the general solution is a lin- 
ear combination of first order Hankel functions of the first and 
second kind, see Eqns. ( |151 1 and (152 1. These mode functions 
are a function of Rk{t), and therefore in the limit of i?^ ^ cx3 
the mode functions approach positive and negative frequency 
modes, while for Rk the modes stop oscillating, and the 
modes exhibit exponentially growing or exponentially decay- 
ing kinematics. 

A simpler way to come to the same answer is to investigate 
Eq. ( 175| l in its limits: These are ujkit) '3> H, and uJk{t) ^ 
H, or equivalently in terms of the frequency ratio: 

R-k(t) ^ 1: It is then possible to write down an approximate 
solutions for Eq. ( 175| l, 
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exp(i /* uJk{t') dt') 



exp(-i /j*^ ujk{t') dt') 



(176) 



(177) 



These are approximately plane waves, but their ampli- 
tude and frequency change as a function of time. This 
ansatz is referred to as the WKB approximation, which 
is valid within the adiabatic limit, when during one os- 
cillation period T — 2n/u)k{t) the relative change in 
the frequency is small (see 1421 ). 
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(179) 



the condition that the ratio Rk(f) is much larger than 
one, which verifies the consistency of adopting the 
adiabatic approximation. 

Hence, the particle production process in the infinite 
past is negligibly small, as all modes oscillate rapidly 
compared to the Hubble frequency. 



R-k(t) <C 1: Here the differential equation ( 175 1, reduces to 

=0. (180) 
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FIG. 4: (Colors online only.) In this figure we compare the quasiparticle production per quantum mode (left column) with its frequency ratio 
(right column), for ts = 1 x 10"^. Parameters are C]VL(t = 0) = 1 x 10^, A^o = 10^ and X = 2 x 10^. The bold plotted dots on the 
left hand side indicate that the frequency ratio is below one, hence the quantum mode corresponds to a super-Hubble horizon mode. While 
on the right hand side we indicated with the bold dots when a change in the mode occupation number is above a certain threshold — here 
ANk 0.004 — to filter out quantum noise fluctuations. 
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The solutions of this equation are exponentially grow- 
ing, ~ exp{Ht), or decaying, ~ exp{—Ht), mode 
functions. 

Therefore, the modes are no longer freely oscillating. 
Instead they get "dragged along" with the spacetime 
fabric. 



Rk(t) = 1: For any mode k there is a time tcrc 

frequency ratio is equal to one, and such that: 



t < 
t > L 



^k{tf > 1 

^k{tY < 1 



when the 

(181) 
(182) 



Therefore it is possible to associate a boundary, the so- 
called ''Hubble horizon", when cj/j = H. For different 
modes with wavenumber k the Hubble horizon occurs 
at different times. Subhorizon modes, i.e., ujk > H sie 
distorted plane wave oscillations, while superhorizon 
modes, i.e., ujk < H, approximately satisfy a harmonic 
oscillator equation with imaginary oscillator frequency, 
and thus they are exponentially growing or decaying. 
The "crossing-time" t = icrossing is referred to as the 
time of "Horizon crossing". The superhorizon modes 
are sometimes said to be frozen modes. Following this 
definition, if we retain the hydrodynamic approxima- 
tion, then in the infinite future all modes will be frozen 
out. 

Please bear in mind that the "Hubble horizon" is dif- 



ferent from the cosmological horizon, see Sec. IV C 3 



One easy way to see this, is to remember that the Hub- 
ble "horizon" is associated with a point in time, which 
is different for each mode, while the cosmological hori- 
zon is the maximum distance signals can propagate for 
any mode. 

The main lesson for the hydrodynamic limit is that it is 
not essential to explicitly solve the differential equation. Of 
course, it is much easier to investigate certain limits of the 
resulting mode functions, instead of investigating the limits 
of the equation of motion. In the hydrodynamic limit both 
methods are accessible. Beyond this limit, where the eikonal 
approximation is valid, we suggest a rather different strategy. 

b. Qualitative behavior of particle production beyond 
the hydrodynamic limit: Ex ante we would like to motivate 
this section by the remark that while the frequency ratio in 
the hydrodynamic limit Rk (t) is a monotonically decreasing 
function in time, the ratio in the eikonal approximation TZk{t) 
is not. Therefore there is some freedom to obtain different 
results to the "conventional" particle production process. We 
demonstrate the correctness of this assertion by referring to 
the numerical simulations reported in HI. 

To obtain a rough estimate on the different qualitative 
regimes of particle production, we use the experience 
gained in the hydrodynamic limit, and simply exchange 



Rk{t) TZk{t)- The eikonal frequency ratio has been 
introduced in Sec. |IV B 3| see Eq. ( |138[ ), as the ratio between 
the modes (modified) frequency uJk (t) — luq \/bk{ t), and the 
(rainbow) Hubble parameter Hk{t), see Eqs. ( 143 i and ( 138 1. 



For early times, when Rk{t) ^ 1 the hydrodynamic and 
eikonal ratios are identical, therefore in both cases approach 
the adiabatic limit, where we expect the particle production 
process to be negligibly small. 

As intimated, the overall slope of the eikonal ratio is not a 
monotonically decreasing function, since 



Uk = 



iuo y/exp{-2Ht) + {k/KY [exp{-2Ht) - 2 (k/K)^] 

-exp{-2Ht) 

(183) 

and the time derivative of the eikonal limit changes its sign at 
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For t < ttuin the slope of the ratio is negative, for t = tturn 
the ratio is given by 
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(185) 



and for t > iturn the ratio is positive. Therefore the eikonal 
ratio has a minimum at tturn, with the maximal particle 
production around this point. After this point the ratio starts 
to increase again, and we shall soon see that the particle 
production process will slow down again. 

To qualitatively describe the particle production process in 
our specific rainbow spacetime, we suggest the following ter- 
minology: 

t — > — oo: At early times almost all modes are "sub-Hubble- 
horizon" modes, and the particle production process 
is negligible. The modes oscillate with much higher 
frequencies as their corresponding Hubble frequencies, 
that is TZkit) > 1. 

t ^ ttuin: As times goes on the mode frequencies are decreas- 
ing, while at the same time the rainbow Hubble fre- 



quencies are decreasing as well, see Sec. IV B 2 Nev 



ertheless, the ratio between them exhibits a minimum at 
^turn, where the particle production process is expected 
to be maximal. Even if the particle production pro- 
cess is maximal, this does not necessarily imply that the 
quantity of particle production is noticeable; the modes 
also need to be "super-Hubble-horizon" modes, or in 
more accurate terminology, we require 7?,fc(tturn) ^ 1- 



t - t, 



r'. If there exists a time t 



tr. 



crossing" ^ liiv.^±v.^ v^^yvuLj ci Liiiiv.^ V — t-crossmgi SUCh that 

^fc(icrossing) ^ 1^ whcrc a modc k crosses the "Hubble 
horizon", there will be a second time t = iic-cntcring, 
where the mode k enters the "Hubble horizon", and 

i?fc(trc-ontoring) ~ 1- We SUggCSt that it is USCful tO 

adopt the following terminology to describe the behav- 
ior of the modes: "freezing of the mode fc" occurs in the 
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\k\ \k\ 

(e) Nf^ (t) projected onto the k-Nf; plane. (f) 7?.fc (i) projected onto the fc-i?fe plane. 



FIG. 5: (Colors online only.) In this figure we compare the quasiparticle production per quantum mode (left column) with its frequency ratio 
(right column), for ts = 1 x 10'^. Parameters are C]VL(t = 0) = 1 x 10^, A^o = 10^ and X = 2 x 10^. The bold plotted dots on the 
left hand side indicate that the frequency ratio is below one, hence the quantum mode corresponds to a super-Hubble horizon mode. While 
on the right hand side we indicated with the bold dots when a change in the mode occupation number is above a certain threshold — here 
ANk 0.004 — to filter out quantum noise fluctuations. 
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time period tcrossing < t < iturn, whereas "melting of 

the mode fc" occurs during iturn < t < trc-cntcring- 

Note, that we have used the approximate relations 

-Rfc(icrossing) ~ 1 and i?fe( tic- entering) ~ 1, instead of the 

more precise i?fc(<crossing) = 1 and i?fc(<ro-entering) = 1- 

This purely technical complication is due to the fact that 
the harmonic oscillator frequency in the eikonal limit, see 
Eq. ( |142| i, can be written as 

nk{t)^ = ujk{tf^Hk{t)^ [l-2exp(2Ft) {k/K)^] , (186) 

and cannot be simplified to ^lk{t)^ — '^k{t)'^ — Hk{t)^- 
Here the exact times for the "crossing" and "re-entering" of 
horizons should be correlated with the sign-change in the 
harmonic oscillator frequency. Nevertheless, it can be seen 
that this does not change the qualitative description for the 
particle process. 



Another novelty in our qualitative understanding of the par- 
ticle production process in our FRW rainbow-spacetime, is the 
connection with the condensed matter point of view: The min- 
imum of the ratio 7?.fe(tturn) for a mode k, see Eqs. ( 183 1, (57 1, 
and (|54li, occurs at 



exp{-2Ht)-2{k/Kf = ^ fc = Wek(t) V2 • 

(187) 

This quantity also appears in the context of conventional con- 
densed matter physics, where it is defined as the crossover 
between the phonon and free-particle region. This borderline, 
the inverse of the heahng, or coherence length ^ li59 J . is given 
by 
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which indicates where each mode starts to decouple from the 
spacetime. In other words, each mode can experience parti- 
cle production, until it becomes free-particle like. Hence, in 
terms of the microscopic physics of a BEC we have a natural 
understanding of 7?.fe(ttuin)- For a detailed treatment of the 
numerical simulations please see |[T|. 

To show the qualitative correlation between the modified 



frequency ratio ( 138 1 with particle production in our specific 
rainbow de Sitter spacetime, we have plotted the ratio for sev- 
eral fc-modes as a function of time, and compared them to 
number occupation plots, see Figs. [4j [5] [6] and [7] This figure 
compares the change in the mode occupation number in each 
mode on the left side, with frequency ratio TZk{t) of the mode 
on the right side, for a number of different scaling times tg. 
Each black line on the left hand side indicates the occupation 
number in the mode A: as a function of time. (The red dashed 
line on the left had side indicates the sudden limit.) On the 
right hand side we have plotted the frequency ratio TZk {t) for 
each of those modes with a different color (online only); grad- 
ually changing from infrared modes (dark red), to ultraviolet 
modes (dark blue). The horizontal (red dashed) line indicates 
where the frequency ratio TZk (t) is equal to one, while the ver- 
tical (blue dashed) line indicates the end of the expansion time 
in our simulations. 



The black dots in the figures to the left indicate when the 
modes cross over from phononic to trans-phononic behav- 
ior, i.e. where each mode starts to decouple from the emergent 
spacetime. We can see that the black dots are located where 



the frequency ratio has its minimum; see Figs. |4(f)] |5(f)[ \6(f) 
|7(f)land[T2(f) 



Figures |4j |5] [6] and [7] correspond to four different scaling 
times ts for ^ 1 x IQ-^, ts = 1 x 10"'', ts = 5 x 10"^ 
and tg = I X 10"^. In all these figures we compare quasi- 
particle production (left column) with the frequency ration 
(right column). The simulation parameters in these simula- 
tions are in correspondence to reference [1|, where the initial 
non-linearity was given by CNL{t = 0) = 1 x 10^, the num- 
ber of atoms in the condensates correspond to A'o = 10^ and 
X = 2 X 10^. 



For example, in Fig. 4(a) we see the quasi-particle produc- 
tion in each mode in the de Sitter region where ts — Ix 10^^. 
The negligibly small population of the modes during, and af- 
ter, the expansion can be explained by means of the ratio plot 
to its right, in Fig. 4(b) As pointed out above, the particle 



production process is large only if TZk <C 1, which is impos- 
sible to achieve for such a relatively large scaling unit tg. We 
also see, that only the first two modes from the bottom of the 
infrared scale cross the "Hubble horizon", such that TZk < 1, 
and (as indicated in the figure) they eventually turn around and 
re-enter the "Hubble horizon" after a few e-foldings. Conse- 
quently, such an experimental set-up is inappropriate for mim- 
icking cosmological particle production. 

In contrast. Fig. |7(a)| shows significant quasi-particle 
production, for the first 100 modes plotted in the figure. Here 
tg — 1 X 10"^, and thus the expansion is two magnitudes 
faster than one we discussed before. The qualitative behavior 
is roughly in agreement with the mode frequency ratio plotted 
to the right in Fig. 



7(b) 



The particle production process 
slows down after each modes crosses over from the phononic 
to free-particle regime (at the blue point), and as we can 
see from the figures, this happens for the four lowest energy 
modes between 5 — 8 e-foldings. 

While the emergent spacetime picture is necessary to un- 
derstand the time-dependent commutator relations — in terms 
of the field operator and its conjugate momentum on a time- 
dependent classical background — the the particle production 
process tends to zero during inflation is naturally explained 
from a condensed matter physics point of view, as it is related 
to the crossover between phononic and trans-phononic quasi- 
particle excitations. 



3. Emergent cosmological horizons? 

The existence of a cosmological horizon in our specific 
emergent spacetime can be investigated by calculating the 
maximum distance, rmax, travelled by small perturbations ini- 
tiated at a certain time t = to and certain point fo. We 
can associate a cosmological horizon to each point (<o,?\)), 
if the maximum distance the signal — here our excitations 
in the Bose-Einstein condensate — can travel is finite. This 
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FIG. 6: (Colors online only.) In this figure we compare the quasiparticle production per quantum mode (left column) with its frequency ratio 
(right column), for ts = 5 x 10"^. Parameters are C]VL(t = 0) = 1 x 10^, A^o = 10^ and X = 2 x 10^. The bold plotted dots on the 
left hand side indicate that the frequency ratio is below one, hence the quantum mode corresponds to a super-Hubble horizon mode. While 
on the right hand side we indicated with the bold dots when a change in the mode occupation number is above a certain threshold — here 
ANk 0.004 — to filter out quantum noise fluctuations. 
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naturally defines a region around the point of emission, and 
its boundary is the cosmological horizon. In the presence of 
cosmological horizons two points in spacetime are said to be 
"causally disconnected" if the distance between them is larger 

than Tniax- 

In "conventional" cosmology cosmological horizons are 
predicted to exist for an infinitely long-lasting de Sitter uni- 
verse, but what is the situation in our spacetime emerging 
from a Bose gas? We would like to address this question 
presently, and show that due to non-perturbative corrections 
there are, strictly speaking, no emergent cosmological hori- 
zons present in the system. 

This can easily be shown, given that one only needs to in- 
tegrate the group velocity of a perturbation emitted at {to, fo) 
to the infinite future. 



lim 

; — ^+00 



dt. 



(189) 



In "conventional" — that is Lorentz-invariant spacetimes — 
group and phase velocities are identical. This is not the case 
for our specific analogue model, where the phase velocity. 
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y^exp{-2Ht) + (k/K)^, (190) 



and the group velocity, 

cgroup :- ^ 



cgexp(-2iJt) + 2{k/KY 
y/clexp{-2Ht) + {k/Kf 



(191) 



are clearly different. 

Within the hydrodynamic limit, when K ^ oo, the two ve- 
locities are equivalent, Cgroup(i) Ck{t). Within this limit all 
signals travel with the same speed, and therefore wave packets 
"keep their shape", i.e. they exhibit no dispersion. In this case 
a signal — sent at {to,fo) — propagating forever will only 
ti-avel a finite distance given by. 



4. Long lasting rainbow inflation 

The present analogue model does not possess horizons in 
the cosmological sense. Consequently we cannot automat- 
ically assume that an infinitely long-lasting de Sitter phase 
would lead to a Planckian spectrum where the temperature is 
connected to the surface gravity at the horizon, see |60|. Thus 
there is no straightforward method to calculate the overall 
temperature for our everlasting de Sitter rainbow spacetime. 

At first, it seems to be a highly contrived question to ask 
for the final spectrum in our emergent rainbow spacetime 
after an infinitely long-lasting inflationary epoch. How- 
ever, we demonstrate by means of numerical evidence — 
from our computer simulations HI — that a sufficiently 
long duration expansion is sufficient to reveal the character- 
istic final shape of the particle spectrum at the end of inflation. 

Before we jump to the result, we would like to show that the 
final Bogoliubov coefficients are indeed time-independent, 
and that the mode occupation number calculated with the 
instantaneous Hamiltonian diagonalization method represents 
"real" quasi-particles, due to well-defined in states in the 
adiabatic regime, and also well-defined out states in the 
free-particle regime. The mode functions in our rainbow 
spacetime interpolate between these two states. 

a. Pair of coupled harmonic oscillators: Given that we 
are only dealing with one region, we can drop the prefactor 
in the equation of motion ( [66| for the auxiliary mode operator 



where 



±k{t) + nk{tfxk{t) = Q, (194) 

o : ifc(t) + (^f™)2xfe(t) = 0, (195) 
o : jck{t) + {ujfYxk{t)^^. (196) 
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In the past cosmological horizons have been repeatedly mis- 
interpreted, and we would like to advise the interested reader 
to carefully read |54|. 

Returning to our problem of investigating the presence of 
cosmological horizons in the Bose gas, the maximum distance 
a signal can travel is given by 



In the infinite past, when the eikonal scale factor approaches 
the hydrodynamic ("conventional") scale factor, the mode 
functions are determined within the adiabatic regime. Thus 
Lo^^^ ^ iOk{t), and the mode functions are approximated by 
Eqs. ( |176| l and ( |177| i. These mode functions represent positive 
and negative frequency modes at every instant of time. More- 
over in the infinite past they represent the adiabatic vacuum 
— for instance see ll28l l42l I6l1 . 

In the infinite future all excitations behave like free parti- 
cles, since 



^max — lim / Cgioup dt. 



OO . 



(193) 



Thus, the emergent spacetime picture strictly does not exhibit 
any cosmological horizons whatsoever. This result can be 
generalized beyond the specific BEC based model, since many 
other analogue models for gravity exhibit super-sonic modifi- 
cations in the dispersion relations will show similar behavior 
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represents the kinetic energy of an object with mass to; 
E = p^/{2m) = h^k"^ /{2rn). (In terms of BEC theory this 
is the "single particle energy".) The time-independence is 
due to the end of the effective expansion in our emergent 
spacetime. The mode functions are the usual positive and 
negative frequency modes, as given in Eq. ( 109| l, where 
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FIG. 7: (Colors online only.) In this figure we compare the quasiparticle production per quantum mode (left column) with its frequency ratio 
(right column), for ts = 1 x 10~^ . Parameters are C]VL(t = 0) = 1 x 10^, A^o = 10^ and X = 2 x 10^. The bold plotted dots on the 
left hand side indicate that the frequency ratio is below one, hence the quantum mode corresponds to a super-Hubble horizon mode. While 
on the right hand side we indicated with the bold dots when a change in the mode occupation number is above a certain threshold — here 
ANk 0.004 — to filter out quantum noise fluctuations. 
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FIG. 8: Schematic description of tlie freezing and melting of quan- 
tum modes. 
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Here we are dealing with well-defined initial and final 
vacua states. This kind of problem can be compared with the 
tanh -expansion, studied for example in ESllMI . In this case 
the Bogoliubov coefficients are — after an infinitely long- 
lasting expansion — time-independent. Please note that this 
analytic calculation does not by itself give any details when 
and how the mode functions change, since the Bogoliubov co- 
efficients for a smooth scale function bk{t) for all t are given 
by the globally defined time-dependent quantities 
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as compared with Eqs. (103 1 and (104i. The Wronskian of 
the mode function is time-independent, therefore it can be 
evaluated at all times. (Note that the general mode functions 
are not necessarily normalized.) Again, Vk and represent 
the in, and and u*^ the out mode functions. (Note that we 
have changed our notation in respect to 1 1 1.) 

By means of the numerical results shown in Figs. (|9|, it ap- 
pears to be rapidly approaching a final asymptotic state. We 
have plotted the final spectra in Fig. 10 so that we easily see 



the relationship between the slope of the line and the inverse 
scaling time t^^. However the temporal duration in our pre- 
vious numerical simulation for tg = \ y. 10^^ has not been 
sufficiently long, as it is not obvious that the particle produc- 
tion process has yet came to an end. In Fig. [7] we can see that 
at the end of the numerical simulation a good fraction of the 
quantum field modes {i.e., roughly |fc| < 40) are frozen. We 
expect particle production in those modes to contribute sig- 
nificantly to the infrared end of the final spectrum. We there- 
fore repeated the numerical simulation for = 1 x lO^'^ (so 
that initially all modes are sub-Hubble horizon modes) with 
two times the previous duration. In addition, the initial non- 



linearity is now C = 2 X 10^ instead of C = 1 x 10^ so 
that all modes are "phononic" at the start of the simulation. 
We also increased our expansion rate from X — 2 x IQ^ to 
10^. As shown in Fig. 11 and in greater detail 



4 X 



X 

in Fig. [121 at the end of the simulation all modes are trans- 
phononic, and the particle production process ceases. From 
our numerical simulations the final particle spectrum does not 
seem to nicely fit a straight line, but it seems conceivable to 
employ standard line-fitting tools to study the final particle 
spectrum as a function of tg and k. (We are currently investi- 
gating this issue.) 



V. CONCLUSIONS AND OUTLOOK 

In this article we put the analogy between a parametri- 
cally excited Bose-Einstein condensate and cosmological 
particle production to the test. Knowing that the analogy 
for mimicking a specific quantum effect in "conventional" 
curved-spacetime quantum-field-theory hinges on the ro- 
bustness of the effect against model-specific deviations, we 
derived the "whole" model-dependent emergent rainbow 
spacetime. Similar work on the acoustic Hawking effect in 
subsonic and supersonic (super)-fluids has been carried out 
in 115J. 



There were two main lessons learnt for the analogue model 
community. First, the specific model we presented — a uni- 
form gas of atoms with time-dependent atomic interactions 
— is in general not robust against the non-perturbative ultra- 
violet corrections. Secondly, we also showed that the anal- 
ogy is sufficiently good for mimicking some aspects of cos- 
mological particle production for finite changes in the size of 
the effective universe. We said "some aspects", because the 
analogy only holds for the low-energy part of the spectrum, 
and therefore the analogy is associated with a certain fc-range. 
These correspond to phononic excitations, bounded by a time- 
dependent parameter |fc| < 1 /^Planck (i)- 

Given that we expect significant deviation from the de- 
sired quantum effect, one might ask the question "Is the 
analogue model we have presented a suitable candidate for 
laboratory experiments?" Previously in |1 1, as well as briefly 
in the current paper, we have presented numerical results 
for cosmological particle production in a "realistic" Bose 
gas. As a matter of fact, despite many possible sources of 
difficulties, for example back-reaction-effects /mode-mixing, 
the phononic regime shows excellent agreement with the 
theoretical predictions |66|. The Bose-Einstein condensate 
enables us to prepare and control a quantum field to such an 
extent, that within a few years time the technology should 
be able to drive "inflation" between two "natural" vacua 
and — that is the outstanding problem |3] |9l |62l — directly 
measure the resulting spectra. Of course, there are also other 
models involving a freely expanding condensate cloud, but 
we suspect similar problems — a growing "Planck-length" 
and the lack of a (strict) cosmological horizon — to 
appear, and would like to stress that those models destroy the 
condensate during the expansion process. We leave it as an 



29 




200 



200 




150 200 



(b) is = 1 X 10" 



1 



-1 

) 

-2 
-3 



-4 








50 



(c) ts=5x 10" 



100 

\k\ 

(d) is = 1 X 10^ 



150 



200 



10 20 30 40 JO 60 70 80 90 100 

I i [ \ I \ I I [ \ \ ] [ ] \ ' I ' l ' ' ' ' ' ' ' 

; I : j I ' ; 1 I ; i i i | i i < i | ■ »— i • i • ■ ■ ■ j • ■ • • 

(e) Duration time T = tji„al - ^initial in percent 

FIG. 9: The time-evolution of the particle spectrum on a single plot, for four different scaling times is. As time goes on the data points 
linearly change their color, changing gradually from violet to turquoise. Note that at low momenta — corresponding to those modes that twice 
"cross the Hubble horizon" by "freezing" and "re-entering" - the final quasi-particle density is approximately log-linear, corresponding to 
rifc ~ Cexp(— (rfc). Some modes never "cross the Hubble horizon", if 7?.fe(iturn) ^ 1. This can be seen in our plots, where we indicated 
^static = [\/23'~^'^*''/\/7qp] ■ ill e,&ch plot with a vertical (red) dashed line. The dashed green line shows the theoretical results for a 
finite de Sitter calculation obtained in the hydrodynamic limit, as presented in Sec. |IV C"T| 




open problem to simulate and check the "robustness" of those 
models. 

So much for the analogy, but what did we learn from the 
deviations occurring in our specific emergent spacetime? Are 
there conclusions to be drawn that are of relevance for the cos- 
mology, or even quantum gravity programme? We leave this 



to our readers, and merely summarize our experience regard- 
ing "trans-Planckian" physics in our emergent spacetime: 

The emergent spacetime we have presented is an exam- 
ple of emergent Lorentz-symmetry. At the infrared end of 
the excitation spectrum we exhibit Lorentz invariance is ex- 
hibited, while non-perturbative corrections from the micro- 
scopic substructure "naturally" break the Lorentz -invariance 
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FIG. 10: The four different particle spectra, corresponding to four 
different scaling times ts, after a finite rainbow de Sitter phase. (Note 
that the "gaps" in the quasi-particle spectrum are a result of the loga- 
rithmic scale that is undefined for (small) negative occupation num- 
bers, a numerical artifact arising in our simulations. For more in- 
formation about the numerical method used to obtain these results 
consult 1 1 1.) 
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FIG. 11: The final particle spectra after a finite duration rainbow de 
Sitter phase for ts = 1 x 10^^ after all modes have decoupled from 
the emergent spacetime geometry. Compare with Fig. |9(d)| 



physics.) The "Planck-length" in our model is correlated with 
the contact potential (scattering length), and thus we are deal- 
ing with a time-dependent "Planck-length". In our expand- 
ing universe the "Planck-length" is growing as well f23], such 
that more and more modes are "trans-Planckian" as time goes 
on. Consequently they gradually decouple from the emergent 
spacetime picture, and their behavior becomes like free parti- 
cles. The rainbow scale factor 

(t) = a(i)/v/l + fc2 

^Planck {ty , (200) 

and the rainbow Hubble parameter 

i/fe(t) = i//(l + fc2£pi^„ckW'), (201) 

are both momentum-dependent. The growth in the "Planck- 
length" forces the rainbow Hubble parameter to approach 
zero, and therefore the universe gradually — mode by mode 
— effectively stops expanding 163] [64l . 

This leads to interesting physics for rainbow inflation: At 
early times all modes are "sub-Hubble-horizon" modes. As 
the expansion goes on some — but not all — modes cross the 
"Hubble-horizon" and become "superhorizon" modes, these 
modes are "frozen modes" and get dragged along with the 
spacetime fabric. Eventually the effective expansion starts to 
slow down — due to the growing "Planck-length" — and a 
process starts that we call "the melting of the modes". After 
a while these modes "re-enter" the "Hubble-horizon" and the 
particle production process is finished. 

However useful these results are to the general relativity 
and cosmology community, we would like to end our conclu- 
sion by commenting on the importance of BEC based ana- 
logue models to the condensed matter physics community. 
There are many aspects, e.g., the time-dependent commuta- 
tion relations, that really seem to require the emergent space- 
time picture to fully understand the physics in our paramet- 
rically excited condensate. The emergent spacetime and the 
Bose-Einstein condensate are two aspects of one and the same 
effect, and one needs to know both points of view to appreci- 
ate the full complexity of this novel state of matter. 



in the ultraviolet regime fTS] [191 |20l El 122 . These correc- 
tions also alter the spacetime picture from "ordinary" space- 
times to the — more "unusual" but conceivable ifTTl — con- 
cept of rainbow spacetimes. These are momentum-dependent 
spacetimes, where the fc-dependence is suppressed in the in- 
frared regime. Further, a time-dependence in the inter-atomic 
potential yields a FRW-type universe for phononic modes, 
and rainbow FRW-type spacetimes for higher-energy modes. 
The borderline that divides the energy scale into phononic 
and trans-phononic intervals, may be viewed as an analogue 
Planck-length. (The physical behavior changes somewhere in 
between the phononic and free-particle regime, in the same 
sense that the Planck-scale is expected to exhibit new laws of 
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(e) (t) projected onto the k-Nf^ plane. (f) 7?.fc {t) projected onto the k-R^ plane. 



FIG. 12: (Colors online only.) In this figure we compare the quasiparticle production per quantum mode (left column) with its frequency ratio 
(right column), for ts = 1 x 10"^. Parameters are CNL{t = 0) = 2 x 10^, A^o = 10^ and X = 4 x 10^. The bold plotted dots on the 
left hand side indicate that the frequency ratio is below one, hence the quantum mode corresponds to a super-Hubble horizon mode. While 
on the right hand side we indicated with the bold dots when a change in the mode occupation number is above a certain threshold — here 
ANk 0.004 — to filter out quantum noise fluctuations. 
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